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ABSTRACT 


The strain pattern analysis (SPA) method is applied to estimate rotor blade 
torsional deflections. The SPA technique requires calculated mode shapes for the tested 
rotor blade and strain measurements from the rotor’s wind tunnel or flight test. The 
Holzer method is developed to calculate the required mode shapes from rotor blade 
stiffness and mass properties and the torsional equation of motion. The Holzer method is 
tested with numerous theoretical and experimental cases and is proven accurate. The 
strain measurements are from wind tunnel tests conducted by the Army, NASA, United 
Technologies Research Center (UTRC) and Sikorsky at DNW with a (1:5.73) model- 
scale UH-60A rotor blade with an advance ratio of 0.301, an advancing tip Mach number 
of 0.8224 and an average Reynolds number of 1,278,729. The SPA method predicts 
slightly larger torsional deflections that compare well with the overall trend and range of 
UTRC static method integrated deflections. The SPA method was evaluated to determine 
the tolerance to change in the number of measurements and the modes applied, errors in 
the measurements, and errors in rotor blade stiffness and mass properties. The method is 


tolerant to all effects except a decrease in the number of measurements and modes. 
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I. INTRODUCTION 


A. ROTOR BLADE TORSIONAL DEFORMATION 


Helicopter rotor blades have high aspect ratios and are highly flexible and 
therefore the effects of aeroelasticity are larger and therefore more important than found 
in more rigid structures. The aeroelastic response of rotor blades includes elastic 
deformation in the flapwise, w, chordwise, v, and feathering, 9, planes of motion. The 
coordinate system for rotor blade deformation in the positive direction is shown below in 
Figure 1. The airfoil shown in Figure 1 is representative of any segment along the rotor 


blade span at rotational speed, Q. [Ref. 1] 





Figure 1. Rotor Blade Coordinate System. 

Aeroelastic deformations of the rotor blade affect helicopter performance, air 
loads, and vibration. Deformation in the feathering plane is known as torsional 
deformation, ¢. The torsional deformation has the most immediate influence on 
aerodynamics since it contributes to the angle of attack, a@, for the rotor blade. The angle 
of attack is the most important factor in the production of lift along the span of the rotor 
blade. Equation 1.1 shows how the inflow angle, ®, and the geometric pitch angle, ©, 
contribute to the angle of attack. These angles are also shown above in Figure 1, with the 


angle of attack shown in gray shading. [Ref. 1] 
a=O0-®@ (1.1) 


Torsional deformation contributes to the angle of attack as a significant term in 
the geometric pitch angle, Equation 1.2, where @, O@ic, and Ois are the collective and 


cyclic pitch control angles, and 6 is the built-in twist. For all terms in Equation 1.2, a 


1 


positive angle denotes the nose up or leading edge up direction in the feathering plane. 
[Ref. 2] 


O(x,t) = Oa(x)+ O(x,t1)+ 00+ 8,. cos @t+6,, sin ot (1.2) 


where 

X = spanwise coordinate 

t = time coordinate 

A subgroup to the geometric pitch angle is the feathering blade angle without 


deformation, , shown in Equation 1.3. The feathering blade angle includes all of the 


terms in Equation 1.2 minus the torsional deformation. [Ref. 2] 


B(x,t) = O2(x)+00+6,, cos at+@,, sin ot (1.3) 


B. MEASUREMENT OF TORSIONAL DEFORMATION 


Given the importance of torsional deformation, a method to measure elastic twist 
during wind tunnel experiments and flight test programs is desired. Unfortunately, 
torsional deformation of a rotating blade is not easily measured during wind tunnel 
experiments or flight-tests [Ref. 3]. The methods to measure or estimate torsional 
deflections during high-speed rotation include the optical projected grid and target 
attitude in real time methods, and the strain-measurement based static and strain pattern 


analysis (SPA) methods. [Ref. 4] 


The projected grid method (PGM) and the target attitude in real time method 
(TART) are both optical methods of measuring deflection. These optical methods are 
only applicable to wind tunnel tests. The projected grid method was developed at the 
German-Dutch wind tunnel (Deutsch-Niederlandischer Windkanal or DNW). With 
PGM, a grid is projected optically onto the surface of the blade and images of the grids 
from the deformed and undeformed blades are compared. Both flapwise and torsional 
deformations can be measured. The method is generally applied over a limited azimuth 
range, y, and extensive video equipment is required. The TART method was developed 
by France’s Office National D’Etudes et Recherches Aerospatiales (ONERA). The 
TART method measures flapwise and torsional deformations for the blade tip at selected 


azimuths. The TART method requires adhesive material embedded with reflective micro 


spheres be applied to the rotor blade. Pulsing light at selected azimuths illuminates the 
micro spheres and the reflected light is recorded by video equipment. Both PGM and 
TART were shown to provide accurate results when compared to the SPA method for 


torsional deformations. [Ref. 4] 


An alternative to optical methods is to use blade strain measurements. 
Detailed strain measurements, if accurate and of sufficient quantity, can be used to 
estimate deformations. One approach, referred to here as the static method, is based on 
simple beam theory and estimates the torsional deformation directly from the measured 
torsion moments, M, and the torsional stiffness term, GJ, as shown in Equation 1.4. [Ref. 


5] 





dg(x) _ M(x) (1.4) 
dx  GJ(x) 


Solving Equation 1.4 for the torsional deformation, @, yields an integral form of 
the static equation shown in Equation 1.5. The non-dimensional radial term, x, increases 
from zero to one at the tip, as shown in Figure 2. The rotor blade extends from the pitch 
link offset, e, to the tip. Integration is done in steps from the tip to the pitch link offset. If 
necessary the torsion moments and torsional stiffness values may be interpolated to span 
locations that coincide with the chosen integral step size. 


_ -M(@) 
A) = | ae 


e 


(1.5) 


pitch link offset=e 


Ax=Ar/R 
te 





x=0 x= 1.0 
Figure 2. Rotor Blade Radial Coordinates. 
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The torsion moment and torsional stiffness values are piecewise uniform terms at 


each individual segment, which provides the descritized Equation 1.6. 


AQ (x) aa (1.6) 


e 


A second approach to estimate deformations from strain measurements is the 
Strain Pattern Analysis (SPA) method, which will be discussed in detail in section C. 
The SPA method uses a least squares estimator derived from fitting calculated torsion 
moment mode shapes to measured torsion moments to estimate the actual deflections 


from calculated deflection mode shapes. [Ref. 6] 


Strain measurements provide an alternative to optical techniques, but have their 
own problems and limitations. To avoid coupling effects, care must be taken to use strain 
gages exclusively for measuring either flapwise or torsional strain [Ref. 7]. The coupling 
between torsional and flapwise deformations occurs when centrifugal forces change blade 
stiffness and mode shapes during rotation [Ref. 7]. When the strain gages successfully 
isolate torsion and flapwise bending mode shapes the results of SPA compared to optical 
methods are shown to be very accurate [Ref. 7]. 


C. STRAIN PATTERN ANALYSIS 


C.J.W. Hassal and D.R. Gaukroger developed the SPA method during the 1970’s 
at the Royal Aircraft Establishment in England. Strain Pattern Analysis was developed to 
confirm mode calculations and to estimate the flap, lag, and torsional deflections of a 
rotating helicopter rotor blade given a set of non-rotating calibration modes and 
corresponding patterns of strain gage measurements from a wind tunnel or flight test. 
The original SPA technique required that calibration modes be determined by loading a 
non-rotating blade and measuring deflections and the corresponding strain gage responses 
[Ref. 6]. William Bousman investigated an extension to the original SPA method 
substituting calculated flapwise mode shapes for the calibration modes [Ref. 3]. This 


thesis applies Bousman’s extension of SPA to the torsional mode shapes. 


The first step in the SPA method is to calculate the least squares estimating 
vector, {q} of size mx/ where m is the number of mode shapes calculated. Equation 1.7 


gives the solution for {q} based on the matrix of calculated torsion moment mode shapes, 
4 


[T], and the vector of strain measurements for one harmonic, {M} [Ref. 8]. The [T] 
matrix is size g x m where g is the number of measurement station locations that must be 


greater than the number of modes calculated, m. The {M} vector is size gx /. [Ref. 6] 


{q} = (TET ETM} (1.7) 
The second step of the SPA method is to confirm the least-squares accuracy of the 


SPA vector. Equation 1.8 solves the least-squares approximation to {M}. [Ref. 8] 


{M} =[T]iqs (1.8) 

The final step of the SPA method is to calculate the estimated deflection vector, 

{@}, from the SPA vector and the matrix of calculated torsional deformation mode 
shapes, [D] as shown in Equation 1.9. The [D] matrix is size n x m, where n is the 
number of calculation stations. The estimated deflection vector is size n x / and is an 


estimate of the actual deflection by harmonic. [Ref. 6] 


{9} =[D]ig} (1.9) 


Both [T] and [D] are obtained analytically, based on blade structural and inertial 
properties. This thesis will focus on developing the application of the SPA method to 
estimate torsional deformation. The accuracy and robustness of the technique will be 
evaluated. The goal of this thesis is to provide the Army/NASA Rotorcraft Division a 
computational program to use in applying SPA to current full-scale and model-scale wind 
tunnel and flight test strain measurements to estimate torsional deformation. 


D. TORSIONAL MODE SHAPE CALCULATIONS 


To apply the SPA method to strain measurements a set of calculated torsion 
moment mode shapes, [T], and a set of calculated torsional deformation mode shapes, [D] 


are required. There are numerous dynamic methods to calculate these mode shapes. 


Modal methods of calculating torsional deformation are accomplished by solving 
the differential equations of motion of the rotor blade [Ref. 9]. The differential equations 
are developed in Chapter II. Exact solutions are available for the differential equations of 
motion representing the non-rotating, uniform blade. For the non-rotating, non-uniform 
or for the rotating case there is no exact solution to the equations of motion. The solution 


can only be approximated, even for the simplest case of a rotating, uniform beam [Ref. 
3) 


10]. The solution obtained with a modal method accounts for the harmonic motion that 
the rotor undergoes in the time domain, as well as the variations that occur spanwise. 
Many methods are available to solve the equation of motion including the Rayleigh 
method, the Rayleigh-Ritz method, the Holzer method, the Myklestad method, and finite 


element analysis. [Ref. 9] 


The Rayleigh method assumes a deflection curve for the first mode. The first 
mode approximation is used to solve the ordinary differential equation for the maximum 
kinetic energy that is set equal to the potential energy. The Rayleigh method achieves 
good approximations for the natural frequencies. The Rayleigh-Ritz method makes the 
first assumed mode a sum of several functions. The functions used to define the first 
mode must collectively satisfy the boundary conditions. The Rayleigh-Ritz method 
provides a more accurate estimate for the natural frequency than the original Rayleigh 


method. [Ref. 9] 


Solution methods known as lumped-mass parameter methods include the Holzer 
method, and the Myklestad method. A lumped-mass model refers to the technique of 
representing a continuous blade as a number of discrete segments connected by massless 
shafts. Each radial location where structural properties are known is considered a 
segment [Ref. 10]. A seven-segment example of a lumped-mass model is shown in 
Figure 3. Each individual segment is located at a discrete span location and can have 
unique values for its mass properties, stiffness properties, and built-in twist angle. A 


uniform lumped-mass model would have identical values for all segments. 





Figure 3. | Lumped-Mass Model of a Cantilevered Beam. 
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The lumped-mass model allows for a partial differential equation of motion to be 
replaced by a set of ordinary differential equations. The Holzer method is used to 
calculate torsional oscillations using a second order differential equation. Myklestad 
extended Holzer’s method to solve for the bending frequencies and modes using a fourth 
order differential equation. The Myklestad method is developed in Professor E. Roberts 
Wood and W. Gerstenberger’s paper titled, “Analysis of Helicopter Aeroelastic 
Characteristics in High-Speed Flight”. [Ref. 9] 


The Holzer method was chosen for calculating the torsional moment mode shapes 
and torsional deformation mode shapes required for this thesis. The Holzer method is 
valid for uncoupled torsion and is directly applicable to my research of uncoupled, 
torsional deformation. A few modifications were required to apply Holzer’s method to 
the rotor blade equations of motion as described in Chapter I]. The Holzer method is 
easily used in tabular form with programs like EXCEL™ when solving for a specific, 
narrowly defined case. A MATLAB™ program was written to compute mode shapes for 
the broadest set of boundary conditions and beam types using the Holzer method. The 
Holzer method is developed more thoroughly in Chapter II. The purpose of this thesis is 
to develop and test SPA for the ability to estimate torsional deformation of a rotor blade 
using calculated, rotating torsional mode shapes and experimentally measured torsion 


moments. 
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Hl. THE HOLZER METHOD 


A. NATURAL FREQUENCY DEFINITION 


A natural frequency occurs where an undamped system vibrates freely and 
continuously with no external force. Equation 2.1 represents the eigenvalue (or 
characteristic value) form of the equation associated with matrices [m] and [k]. [Ref. 5] 

[A]{u} = o° [m]{u} (2.1) 

Equation 2.1 possesses a nontrivial solution when the determinant of the 

coefficients of u; vanish. The nontrivial solution is expressed in Equation 2.2 known as 


the characteristic or frequency equation. [Ref. 5] 


A(@’) = 





2 
k,-@ m,, 





=0 (2.2) 
where 
A(q@”) is the characteristic determinant 


The characteristic equation is of degree p in w’, and possesses in general p 
distinct roots that are called characteristic roots or eigenvalues. The p roots are denoted 
Q, QO: ,...0, and the square roots of these quantities are the system natural 
frequencies @, (m=1,2,...p). The natural frequencies can be arranged in order of 
increasing magnitude where @, <@, <...<@,. The lowest frequency @, is referred to as 
the fundamental frequency, and for many practical problems it is the most important 


natural frequency. [Ref. 5] 


Associated with every one of the natural frequencies @,, is a certain nontrivial 


vector {uv} (m=1,2,...p) of real numbers uw,,, where {u} is a solution of the eigenvalue 


such that Equation 2.3 is valid. [Ref. 5] 
[k]{u},, =n Lme]{u},, (2.3) 


where 


m=1,2,....p 


The vectors {wu} are known as characteristic vectors or eigenvectors. These 


vectors are unique only in the sense that the ratio between any two elements w,,, and w,,, 


is constant. The value of the elements themselves is arbitrary. We say that the shape of 
the natural modes is unique, but the amplitude is not. [Ref. 5] 


B. THE ORIGINAL HOLZER METHOD 


For the Holzer’s method, calculations are initiated by assuming a value for the 
frequency and applying a unit deformation to the blade tip with zero value for tip 
torsional moment. Then proceeding from tip to root, the torsion moments and torsional 


deformations of each segment are calculated. [Ref. 5] 


The partial differential equation (PDE) of motion for torsional inertia is shown in 


Equation 2.4 for a simple, non-rotating blade modeled as a non-uniform beam. [Ref. 5] 


OM (x,t) 0° (x,t) 
———= = I(x) 2.4 
yy (x) ve (2.4) 
The torsional mass moment of inertia term, I, is a function of mass and geometry. 
The Torsion Moment, M, in Equation 2.4 is replaced with the time domain equivalent of 
Equation 1.4. For the lumped-mass model the torsional stiffness and the mass moment of 
inertia terms are piecewise uniform at each segment. The result of applying the above 


terms to Equation 2.4 is seen in Equation 2.5. [Ref. 5] 


- o(x.1) <9 Tae t) (2.5) 
ox ot 





For free vibration the motion is simple harmonic and substitution for the torsional 
deformation term can be done to separate the time and space variables. Separation of 
variables for torsional deformation with simple harmonic motion of frequency, @, is 


shown in Equation 2.6. [Ref. 1] 
P(x,t) = e'P(x) (2.6) 


The variable, e, in Equation 2.3 is the base-e exponential and is not to be confused 
with the pitch link offset term. Substituting for the torsional deformation in Equation 2.5 


with Equation 2.6 and carrying out the differentiating and canceling out the imaginary 


10 


terms we solve for the ordinary differential equation (ODE) shown in Equation 2.7. [Ref. 


5] 


=-1@°9(x) (2.7) 





GJ £9) 
ox 


The right side of Equation 2.7 represents the torsional inertia at any segment. The 
negative sign in front of the torsional inertial term represents a negative direction. In 
Figure 4 the torsional inertia term is shown in the negative direction according to the 


right hand rule 


An analysis of a single lumped-mass segment develops the incremental solution 
to the total torsion moment at any segment. The free body diagram for the single 


segment in Figure 4 helps derive the total torsion moment solution. [Ref. 11] 


i T atl 
+ 
lw} 
Figure 4. Free Body Diagram of a Lumped-Mass Segment. 


Figure 4 shows the total torsion moment of a segment will equal the segment’s 
torsional inertia plus the torsion moment from the right of the segment as represented in 


Equation 2.8. The segment number is represented by the variable, n. [Ref. 9] 


eae 


n+l 


+1,0°@, (2.8) 
Equation 2.9 represents the torsional deformation at a single segment and 


Equation 1.6 is substituted for the term, Ag,,,. The term Ax,,, in Equation 2.9 is equal 


to the difference between the span location of segments, n+1 and n. [Ref. 5] 


Tr, AX 


?, = Prat AQ, 41 = Pra = ea (2.9) 


n+l 


Equations 2.8 and 2.9 combined with the boundary conditions for the lumped-mass 


system establishes the solution equations for Holzer’s method. 
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The presented example of the Holzer method is a twenty segment lumped-mass 
system. The first step to obtain a solution using the Holzer method is to set the torsional 
deformation at the last segment equal to one. Segments will always be numbered from 
one at the root and will ascend in order to the last segment for all lumped-mass model 
applications in this thesis. Equation 2.10 represents the first step of the Holzer method. 
[Ref. 12] 


¢,, = 1 radian (2.10) 


The second step of the Holzer method is to solve for the total torsion moment at 
the last segment by substituting zero for the free end torsion moment and one radian for 
the last segment’s torsional deformation into Equation 2.8. Equation 2.11 represents the 


result of the second step of the Holzer method. [Ref. 12] 
Ty = In” (2.11) 


The third step of the Holzer method is to solve for the next to last segment’s 
torsional deformation using Equation 2.9. The result of step three for a twenty-segment 
system is shown in Equation 2.12. [Ref. 12] 


Ty Ary 
GIa3 





g,, =1- (2.12) 


The fourth step of the Holzer method is to solve for the total torsion moment at 
the next to the last segment using Equation 2.8. The result of step four is shown in 


Equation 2.13. [Ref. 12] 
Tig = Ty) + 1150 0 (2.13) 


The fifth step of the Holzer method requires that steps three and four be repeated 
for each segment in descending order until the solution for segment one at the root is 
obtained. It is important to note that the entire iteration from tip to root is done once for 
each test frequency. The sixth and final step of the Holzer method is to check the root 
solution against the root boundary conditions. The natural frequencies exist where the 


root solution and root boundary conditions match. [Ref. 12] 
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The set of each segment’s torsion moment solution can be plotted versus span to 
give the torsion moment mode shape. Each natural frequency will have an independent 


set of mode shapes. [Ref. 5] 


Three possible root boundary conditions are used for the applications described in 
this thesis. The first condition is for a Free-Free beam. If no external torsion moment is 
applied the torsion moments at both the right and left ends of the beam are zero. 
Additionally, the torsional deformation at the right end is set at one radian. The second 
condition is for a Fixed-Free cantilevered beam where the left end of the beam is clamped 
and therefore has a torsional deformation equal to zero. The right end of the beam is free 
and has a torsion moment of zero and a torsional deformation set at one radian. The third 
condition will be called Flex-Free, where the left end is supported by a flexible root 
restraint that is represented as a spring of stiffness, kroot. The right end is free and has a 
moment of zero and a torsional deformation set at one radian. [Ref.5] The left end 


boundary condition is represented by Equation 2.14 [Ref. 10]. 


Proot = K poo € p / GJ (2. 14) 


‘oot 


The natural frequency solutions can be seen graphically by plotting the root 
boundary condition term of torsion moment or torsional deformation against the range of 
test frequencies. The natural frequencies for the Free-Free case are where the root torsion 
moment equals zero. The natural frequencies for the Fixed-Free case without a flexible 
root are where the root torsional deformation equals zero. The natural frequencies for the 
Fixed-Free case with a flexible root are where the root torsional deformation equals a 


positive value of, ¢,.,,, for odd numbered modes (i.e., the first, third, and fifth natural 
frequencies) or a negative value of, @,,,,. for even numbered modes (1.e., the second, 
fourth, and sixth natural frequencies). [Ref. 2] 

The Holzer method can be represented in transfer matrix form as shown in 


Equation 2.15. To utilize the transfer matrix method the torsional deformation and total 


torsion moment at the last segment must be solved for prior to using the matrix. [Ref. 2] 


(| ; AX, / CS 4 
sa ; (2.15) 
Ty, (1,0 (AX CTn VO) AW ID et 


nt+1 
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C, ROTOR BLADE TORSIONAL EQUATION OF MOTION 


To obtain accurate natural frequencies and mode shapes with Holzer’s method, 
the torsional equation of motion for the system must be specific for rotating blades. The 
development of coupled, non-uniform, rotating blade equations of motion for flapwise, 
chordwise and torsion are based on the work of Hodges and Dowell [Ref. 13]. The 
ordering scheme eliminated all of the third order terms except terms involving rotation 
and torsional deformation. To obtain the uncoupled torsion equation the applied moment 
equals zero for the modal application, the flapwise and chordwise deformations were set 
at zero, the tension term and torsional mass moment of inertia terms were substituted 
according to Hodges and Dowell, and third order terms were eliminated. The torsional 


partial differential equation of motion is shown in Equation 2.16. [Ref. 2,13] 


2 2 
(GJ + Ch?) OD nk! OED inky ~k,,?}Q?6(x,t)cos2B=0 (2.16) 
IX 


where 


km = the polar moment of inertial, k,, =, Kn, + Kn, 


= _ _ 1 (mass) (27 p2 _ 2 
C= Tension = 2 (Ar) (Q°)(R° -r-) 


ka = Radius of Gyration of blade cross sectional area 





The blade cross-section thickness is much less than the chord. Thus, the mass radii of 


gyration simplify as follows in Equation 2.17. [Ref. 2] 
Kg, 2 Km > Kg * Ky (2.17) 


m, 


The torsional mass moment of inertia, I, is equal to mk,” and yields Equation 


218. 





0° (x,t) _ 2 Gxt) 
ax? or? 


The tension term, radius of gyration, torsional stiffness, torsional mass moment of 


(GJ +Ck,) ~1279(x,t)cos2B =0 (2.18) 


inertia and feathering blade angle are all piecewise uniform at each segment of the blade 


when it is modeled as a lumped-mass system. The first term of Equation 2.18, 


GJ(0°¢/dx’), is the same as the torsional stiffness term in Equation 2.5. The second 


term of Equation 2.18, Ck 1 (a°¢/ dx”) is called the tension-torsion term and tends to 


14 


untwist a pre-twisted blade due to centrifugal force. The third term in Equation 2.18, 
f (0°¢/ dt’), is the same as the torsional inertia term in Equation 2.5. The fourth term in 
Equation 2.18, JQ’¢cos2Z, is called the “tennis racket” term and also tends to untwist a 
pre-twisted blade. The untwisting effect from the tension-torsion term and the “tennis 
racquet” term stiffens the rotating blade as a function of Q*. Separation of variables is 
done by substituting for the torsional deformation in Equation 2.18 with Equation 2.6. 
Differentiating @ and canceling the imaginary terms produces the simplified equation. 
Equation 2.19 shows the resulting torsional ordinary differential equation (ODE) of 


motion for rotating blades. [Ref. 13] 





(GJ+Ck,’) 7H) + lo’ @(x) — 12’¢(x)cos2B = 0 (2.19) 
IX 


0 


Equation 2.19 is rearranged to facilitate use with Holzer’s method. Equation 2.20 
gives the Holzer form of the rotating blade equation that will be applied in this thesis. 
The tension-torsion term and the “tennis racket” term are the only differences when 


compared to Equation 2.7. 


2 
(GJ + ck, SE = —I¢(@" — Q? cos2f) (2.20) 
IX 
D. |THE HOLZER METHOD APPLIED TO ROTOR BLADES 


The following modifications are required to apply Holzer’s original method using 
the rotor blade torsional equation of motion shown in Equation 2.20. In place of 
Equation 2.8 we apply the same theory to the free body diagram in Figure 5 with the 


added “tennis racket” term from Equation 2.20. The result of this analysis is shown in 


Equation 2.21. 
16% 2 
T. dcos2p 
— —> T 
—> n+l 
lw 
Figure 5. Free Body Diagram of a Lumped-Mass Rotor Segment. 
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, = Thuy + 1,0,(@" -Q? cos2B,) (2.21) 


Equation 2.22 accounts for the tension-torsion term in the dynamic analysis of 


rotor blades and is used in place of Equation 1.6. 


TAx 


Aé = ——_____ 
? (GJ +Ck,’) 


(2.22) 


Equation 2.23 replaces Equation 2.9 and represents the torsional deformation at a 


single segment where Equation 2.22 is substituted for the term, A@,,,. 


TAG 


(GIyu, + (Ch nat) 





Pn = Pus —APrst = Pra — (2.23) 

Equations 2.21 and 2.23 combined with the boundary conditions for the rotor 
allow a solution using Holzer’s method. The first step of the rotating blade Holzer 
method is the same as in the original method with the torsional deformation at the last 
segment equal to one radian as shown in Equation 2.10. The second step of the Holzer 
method solves for the total torsion moment at the last segment by substituting the 
boundary condition, T2; = 0, and Equation 2.10 into Equation 2.21. The result of step 
two is shown in Equation 2.24. If a twenty segment lumped-mass model of a rotor blade 


is used, then. 
Ty) = 1) (@ —Q’ cos2B,y) (2.24) 


The third step of the Holzer method for the rotating blade solves for the next to 
the last segment’s torsional deformation using Equation 2.23. The result of step three is 
shown in Equation 2.25. 


= I, AX 
(Gy) +(Cky a9) 





Ay =1 (2.25) 


The fourth step of the Holzer method for a rotating blade solves Equation 2.21 for 
the total torsion moment at the next to last segment. The result of step four is shown in 


Equation 2.26. 


Lig = bg 1,0,)(@° —Q’ cos 2B») (2.26) 
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The fifth step of the Holzer method for rotating blades repeats steps three and four 
for each segment in descending order until the root solution is obtained. The sixth step 
for rotating blades is the same as in the original method and requires that the root solution 


be compared against the root boundary conditions. 


The Holzer method for rotating blades can be represented in transfer matrix form 


as shown in Equation 2.27. [Ref. 2] 


par t (Chg nat) 


_ Ax,.1,(@° — Q’ cos2B,) BA (| (2.27) 
I,(@° —Q? cos2Z,) (GJ, +(Ck,) ne 


n+l 


1 ~ Ax, (GI 





n+l 


Steps one through six of the Holzer method as presented for rotor blades is written 


into a MATLAB™ routine listed in Appendix A. 
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il. EVALUATION OF THE HOLZER PROGRAM 


A. INTRODUCTION 


The Holzer program developed here is validated using a sequence of increasingly 
more complex test cases. These test cases include a mixture of theoretical and 
experimental tests. The theoretical test cases start with the simplest case of a uniform 
beam and increase in complexity until a non-uniform, rotating beam is evaluated. 
Rotating and non-rotating evaluations are done for both the Fixed-Free and Free-Free 
boundary conditions to demonstrate that the program is valid for these applications. The 


theoretical test cases will be presented in the following order: 

e Test Case la: Uniform, Fixed-Free, Non-Rotating Beam 

e Test Case 1b: Uniform, Fixed-Free, Rotating Beam 

e Test Case 2a: Free-Free, Non-Rotating, Uniform Beam 

e Test Case 2b: Free-Free, Non-Rotating, Non-Uniform Beam 
e Test Case 3a: Non-Uniform, Fixed-Free, Non-Rotating Wing 
e Test Case 3b: Non-Uniform, Fixed-Free, Rotating Wing 


The accuracy of the Holzer program is also demonstrated using a set of realistic, 


experimental cases. The experimental test cases will be presented in the following order. 
e Test Case 1: Non-Uniform, Non-Rotating, Free-Free, Full-Scale Blade 

e Test Case 2a: Non-Uniform, Non-Rotating, Flex-Free, Model-Scale Blade 

e Test Case 2b: Non-Uniform, Rotating, Flex-Free, Model-Scale Blade 


All of the theoretical and experimental test cases are solved with the Holzer 
program using Equations 2.21 and 2.23. For the non-rotating cases the equations reduce 


to the form of the original Holzer method’s Equations 2.8 and 2.9. 
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B. THEORETICAL TEST CASES 
1. Uniform, Fixed-Free Test Cases 


a. Non-Rotating 


For the simple case of a uniform, cantilevered beam the Fixed-Free 
boundary conditions apply. This case can be solved exactly. The approximate solution 
from the Holzer method will be compared with the exact solution to validate the Holzer 
code. The Holzer method will be applied to four different, equivalent lumped-mass 
models of a six-foot long beam. Each model will use a different number of lumped-mass 
segments to analyze the effects of discretization on the Holzer solution. Equation 3.1 will 
be used to solve for the exact solution of the non-rotating, uniform rotor blade natural 
frequency. Equation 3.2 is the equation for the exact solution of the torsional 


deformation mode shape. [Ref. 1] 


1 GJ 
Oop = (m——)z,|— ol 
NR ( > IR? ( ) 
where 
,, = the non-rotating natural frequency 
m = the mode number 
m m+ : 1 
Orie = (-1"") sin((m — a (3.2) 


where 


Pvp = the non-rotating torsional deformation mode shape 

The uniform beam tested for this case is six feet long, has a torsional mass 
moment of inertia of 40 Ib-ft-sec* and a torsional stiffness of 3x10° Ib-ft?. To model a 
uniform beam into an equivalent system of n evenly spaced segments the total length of 
the beam is divided by (n-1) to get the span length between each segment, Ax,. The 
equivalent segment torsional mass moment of inertia equals the original torsional mass 
moment of inertia multiplied by Av,. The torsional stiffness, GJ, remains unchanged in 
the conversion to an equivalent lumped-mass system. Table 1 shows the properties for 
four cases using different numbers of elements. For each lumped-mass model (i.e., 
seven-segment model, 13-segment model, 25-segment model, and 49-segment model) the 
properties listed are identical for each segment of the same model. 
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—— Models “ "n" Segments 


GJ 
lb-ft-sec” lb-ft? 


Table 1. Structural Properties for Uniform, Non-Rotating, Fixed-Free Beam. 








rota, 
=e 
| 25 
fe ei 


The Holzer program was run four separate times, once for each lumped- 
mass model of the uniform beam and the first three natural frequencies and mode shapes 
were computed. The calculated natural frequencies are listed in Table 2. Equation 3.1 
was used to calculate the exact solution for the first three natural frequencies and is the 


basis for the error values listed in Table 2. 


Model error error error 
a7 sec a sec a. sec 


Se 


Pexact | m7 [pasa 8s | 





Table 2. Natural Frequency Results for the Uniform, Non-Rotating, Fixed-Free Beam. 
As expected, the error decreases as the number of segments increase. 
Errors are inherent with approximation solutions such as the Holzer method. The errors 
are greater for the higher natural frequencies. An error of one percent is achieved for the 


49-segment case, for the first natural frequency solution. 


These calculations are illustrated in Figure 6 that shows the first three 
exact and Holzer calculated natural frequencies versus the number of segments used. 
Each exact solution is shown as a horizontal line at its natural frequency. Each natural 
frequency also has four different attempts to calculate the natural frequency using models 
with an increasing number of segments. Each calculated natural frequency is marked 
with a diamond symbol and the number of segments. The lines between calculated 
solutions simply connect the data points. For all three natural frequencies the calculated 
values converge towards the exact solution with an increasing number of segments. 
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Third Mode = 358.5 rad/sec 


Second Mode = 215.1 rad/sec 


Natural Frequency (rad/sec) 


First Mode = 71.7 rad/sec 





Oo 5 10 15 20 25 30 35 40 45 50 
Segments 


Figure 6. | Natural Frequency Accuracy by Segment Number. 

Figure 7 displays a plot of the Holzer program natural frequency solutions 
for the 49-segment model. The y-axis is the root (segment one) torsional deformation. 
The x-axis is the range of test frequencies attempted. The natural frequencies are where 
the root boundary condition of zero torsional deformation is met. The natural frequency 
solution data points are marked with a square. Each root torsional deformation solution 
is plotted against the respective test frequency at an increment of five (rad/sec). The 
plotted line connects each data point for test frequencies ranging from zero to 400 


(rad/sec) as chosen during program execution. 
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O root B.C. met 
— full set of root solutions 


=0 


Torsinal Deformation @ Root (rad) rR 





Oo 50 100 150 200 250 300 350 400 
Frequency (rad/sec) 


Figure 7. Natural Frequency Graphical Solution to the 49-segment, Uniform, Fixed-Free, 
Non-Rotating Beam. 


Figures 8, 9, and 10 plot the modes shapes for the first, second, and third 
natural frequencies respectively for the 49-segment Holzer solution. The Holzer 
calculated torsional deformation is plotted against span location. Each of the 49- 
segments has its own torsional deformation value marked by a square. The solid line is 


the exact solution plotted using Equation 3.2. 
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Torsional Deformation (rad) 


Torsional Deformation (rad) 





“0 oA 


Figure 8. 





“O O41 


Figure 9. 


Natural Frequency = 71.0 rad/sec 


O calculated mode shape 


0.2 0.3 o4 o5 0.6 O7 08 og 
oR 


49-Segment First Mode Shape. 


Natural Frequency = 212.8 rad/sec 


O calculated mode shape 


0.2 0.3 o4 os 0.6 O7 08 og 
oR 


49-Segment Second Mode Shape. 
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Natural Frequency = 354.4 rad/sec 


O calculated mode shape 


Torsional Deformation (rad) 





oR 


Figure 10. 49-Segment Third Mode Shape. 

All three mode shape plots show excellent agreement with the exact 
solution. Figure 11 shows the deviation of the 49-segment calculated mode shapes with 
the exact solution mode shapes versus span for the first three modes. Each data point is a 
single deviation between the exact and calculated torsional deformation by span location. 
The lines simply connect the data points. The largest error is observed in the third mode 


on the center portion of the beam, approximately -0.065 (tip deflection is 1.0). 
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| eel 
—*- First Mode 
—& Second Mode 
—& Third Mode 


Deviation (rad) 





in) o4 0.2 0.3 o4 05 0.6 O7 08 og 4 
oR 


Figure 11. | Mode Shape Deviation of the 49-Segment, Calculated versus Exact. 
Based on the results presented the Holzer program is valid for the uniform, 


Fixed-Free, non-rotating test case. 


b. Rotating 


The exact solution for the Fixed-Free, uniform, rotating natural frequency 


is represented by Equation 3.3 [Ref. 10]. 


Dror = V (yp. t Q’) (3.3) 


This test case is run with the same structural properties used in the 
previous non-rotating test case listed in Table 1. The rotational speed is 150 (rad/sec) and 


the natural frequency results from the Holzer program are shown in Table 3 


Model error error error 
ay sec ay sec at sec 


fe 


i3_[ 165.0 | 08 | 2542 | 31 | 3703 | 47 
P Exact [1663 | | 202.2 | | 3886 |_| 





Table 3. Natural Frequency Results for the Uniform, Rotating, Fixed-Free Beam. 
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The exact rotating natural frequency solutions in Table 3 are found by 
using the exact non-rotating natural frequency solutions from Table 2 with Equation 3.3. 
As with the non-rotating case, the most accurate solutions are for the first mode, with the 
errors increasing with increased mode number. The third mode is still accurate to one 
percent. It is unknown why the rotating modes are more accurate than their non-rotating 


counterparts. 


The natural frequency error decreases with an increase in segments used. 
Figure 12 shows the convergence of the Holzer natural frequencies towards their 
respective exact solutions as the number of segments increase. Each natural frequency 
solution is marked with a diamond and the number of segments used in the calculation. 


The lines connect the data points. 


Third Mode = 386.6 rad/sec 


Second Mode = 262.2 rad/sec 


Natural Frequency (rad/sec) 





Segments 


Figure 12. Rotating Natural Frequency Accuracy by Segment Number. 
Rotating torsional mode shapes are theoretically identical to the non- 
rotating mode shapes [Ref. 2]. In reality there is some variance between rotating and 
non-rotating mode shapes. Figure 13 plots the deviation between the rotating and non- 


rotating mode shapes versus span for the first three modes. Each data point is the single 
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point deviation for an individual segment. The lines connect the data points for each 


mode. 


Deviation (rad) 


—* First Mode 
— 6 Second Mode 
—6 Third Mode 





“o o4 0.2 0.3 o4 O05 0.6 O7 0.8 og 1 
oR 


Figure 13. Mode Shape Deviation of the 49 Segment, Rotating versus Non-Rotating. 
The deviation for the first mode shape is large in comparison with the 
second and third mode shape deviations. Despite the relatively large deviation for the 
first mode, the scale is in the thousandths and the largest deviation is less than -0.003 (tip 
deflection is 1.0). The Figure 13 mode shape deviations are between Holzer calculated 
rotating mode shapes and Holzer calculated non-rotating mode shapes. The larger 
deviations presented in Figure 11 are between the Holzer calculated non-rotating mode 
shapes and the exact mode shapes. Based on the results presented the Holzer program is 
valid for the uniform, Fixed-Free, rotating test case. 
2, Free-Free, Non-Rotating Test Cases 


a. Uniform 


A three-segment model of a uniform, Free-Free beam with structural 
properties listed by segment in Table 4 is tested with the Holzer program. This test case 


has no exact solution or published solution to compare the calculated results with. This is 
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a reduced or simplified case from Thomson’s non-uniform, Free-Free example 12.4.1 that 


is discussed in the next section. [Ref. 11] 


3-Segment Lumped-Mass Model 


Segment I Ar GJ 
Number kg-m* m N-m? 


Pt | 22 | 0 | 8.0E+06] 


Table 4. Structural Properties for Uniform, Non-Rotating, Free-Free Beam. 








Figures 14 and 15 plot the mode shapes for the first two modes. The data 
points represent the torsional deformation at each of the three segments. The line simply 


connects the data points. 


Natural Frequency = 95.3 rad/sec 


—-6- calculated mode shape 


Torsional Deformation (rad) 
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Figure 14. Free-Free, Uniform First Mode Shape. 
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Natural Frequency = 165.1 rad/sec 


—-6- calculated mode shape 


Torsional Deformation (rad) 


in) O41 0.2 0.3 o4 as 0.6 O7 
oR 





0.8 og 1 


Figure 15. Free-Free, Uniform Second Mode Shape. 


The results of this uniform reduced case from Thomson are consistent 


with the non-uniform case presented next. The Holzer program is valid for the uniform, 


Free-Free case. 


b. Non-Uniform 


This test case is Thomson's example 12.4.1 for a slightly non-uniform, 


Free-Free beam. Table 5 shows the structural properties for the individual segments. 


[Ref. 11] 


3-Segment Lumped-Mass Model 


Segment I Ar GJ 
Number kg-m* m N-m 


aa 


Table 5. Structural Properties for Non-Uniform, Non-Rotating 





ea es 06 
2 





, Free-Free Beam. 


The natural frequency solutions are shown in Figure 16 where the root 


boundary condition for the Free-Free case is a torsion moment equal to zero. The natural 


frequencies are marked with a square. The line simply connects each test frequency's 
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root torsion moment result at test frequency increments of five (rad/sec) from zero to 300 


(rad/sec). 


2000 





OQ root B.C. met 
— full set of root solutions 


= 
pail 
o 
oa 


1000 


Total Torsion Moment (N-m) at segment 1 (*10-3) 


in) 50 100 150 200 250 300 
Frequency (rad/sec) 


Figure 16. | Natural Frequency Graphical Solution to the Free-Free, Non-Uniform Beam. 
The Holzer program gives the exact answers found in the Thomson 
textbook for frequencies and mode shapes. The first natural frequency is at 123.6 
(rad/sec) and the second natural frequency is at 202.5 (rad/sec). Figures 17 and 18 show 
the first and second mode shapes. The squares mark each segment’s deflection and the 


line connects the data points. 


Bal 


Natural Frequency = 123.6 rad/sec 


—6 calculated mode shape 


Torsional Deformation (rad) 
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Figure 17. |. Non-Uniform, Free-Free First Mode Shape. 


Natural Frequency = 202.5 rad/sec 


Torsional Deformation (rad) 





—-6- calculated mode shape 
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rR 


Figure 18. Non-Uniform, Free-Free Second Mode Shape. 
This test case is only slightly non-uniform and therefore the mode shapes 


shown in the previous test case can be compared to Figures 17 and 18 for a general 
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comparison. When Figures 14 and 15 are compared to Figures 17 and 18 the mode 
shapes follow the same general trend. Based on the presented results the Holzer program 
is valid for the non-uniform, Free-Free case. 

3. Non-Uniform, Fixed-Free Test Cases 


a. Non-Rotating Wing 


This test case is Scanlan's example 7.4. The structural properties are listed 
in Table 6 are for each segment. This example is a model of a fixed-wing aircraft. [Ref. 
12] 

7-Segment Lumped-Mass Model 
Segment I Ar GJ 
Number | lb-ft-sec” ft lb-ft” 
Po | | 4.0 +09 
4.0E+09 


Peas) 40 | E09 


Table 6. Structural Properties for Non-Uniform, Non-Rotating, Fixed-Free Wing. 








The Holzer program gives the exact answers found in the Scanlan 


textbook for natural frequencies and mode shapes. The natural frequency solutions are 


listed in Table 7. 
@ 1) 03 
rad/sec rad/sec rad/sec 


Table 7. Natural Frequency Results for the Non-Uniform, Non-Rotating, Fixed-Free Wing. 









Figures 19, 20 and 21 plot the first three mode shapes. The individual 
segment deformation value is marked with a square. The line simply connects the data 


points. 
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Natural Frequency = 98.8 rad/sec 
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Figure 19. Non-Rotating Wing First Mode Shape. 


Natural Frequency = 351.9 rad/sec 


-6- calculated mode shape 
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Figure 20. Non-Rotating Wing Second Mode Shape. 
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Natural Frequency = 439.1 rad/sec 


—6 calculated mode shape 


Torsional Deformation (rad) 
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Figure 21. _Non-Rotating Wing Third Mode Shape. 
The Holzer results for the non-rotating, non-uniform, Fixed-Free case of 
the Scanlan wing are exact with the published solution. The Holzer program is valid for 


this test case. 


b. Rotating Wing 


This test case takes the structural properties listed in Table 4 and applies a 
rotational speed of 150 (rad/sec). The rotating natural frequencies from the Holzer 
program are identical to the exact solution. The exact rotating natural frequency solution 
is obtained using Equation 3.3 and the non-rotating natural frequencies from Table 7. 


Table 8 lists the rotating natural frequency results. 


@ 1) 003 
rad/sec rad/sec rad/sec 


Table 8. | Natural Frequency Results for the Non-Uniform, Rotating, Fixed-Free Wing. 









Theoretically the mode shapes in rotation are identical to those for the 
non-rotating condition. In reality there is some deviation between rotating and non- 


rotating mode shapes. Figure 22 presents the mode shape deviation between the rotating 
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and non-rotating mode shapes for the first three modes. Each symbol marks the deviation 


at an individual segment and mode. The lines connect the data points for each mode. 


ceaiaieeienieeeal 
—*- First Mode 
—@ Second Mode 
—-& Third Mode 


Deviation (rad) 





oR 


Figure 22. Deviation of Rotating versus Non-Rotating Wing Mode Shapes. 

The deviations in Figure 22 are extremely small with the largest value 
being less than 0.0004 radians. Based on the presented results the Holzer program is 
valid for the non-uniform, Fixed-Free, rotating test case. 

C. NON-UNIFORM ROTOR BLADE TEST CASES 
1. Full Scale UH-60 Rotor Blade 


a. Overview of Experiment 


In November 1990, Engineers at the NASA-Ames Research Center 
measured dynamic properties of a full-scale production and instrumented UH-60A Black 
Hawk rotor blades. NASA measurements from the shake tests of a production blade will 


be compared against the Holzer program results. [Ref. 15] 


The UH-60A rotor blade is 24.15 feet in radius and is composed of a 
titanium spar, a honeycomb core, fiberglass skin, and an aluminum tip cap. The root 


offset is 2.65 feet. [Ref. 15] 
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b. Non-Rotating, Free-Free Natural Frequencies and Mode Shapes 


Table 9 shows the structural properties for the full-scale UH-60A rotor 
blade from a NASA database [Ref. 16]. The list of structural properties does not include 
terms required to generate the tension-torsion or the “tennis racket” terms because these 


terms are only a factor for rotating test cases. 


af 


GJ | r GJ | r 
lb-ft 7 lb-ft-sec 7 ft lb-ft 7 lb-ft-sec 7 ft 


Table 9. | _76-Segment Full-Scale UH-60A Blade Structural Properties. [After: Ref. 16] 
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Table 10 presents the natural frequency results and errors in comparison to 


the NASA measurements. 









| NASAMeasurement | 46.6 | | 85.1 | 


Table 10. Natural Frequencies for the Full-Scale UH-60A Rotor Blade. [After: Ref. 15] 
The result for the first natural frequency barely exceeds one-percent 
compared to the experimental result. The error for the uncoupled Holzer calculated 
second natural frequency is significant at greater than nine-percent. The results are 
acceptable since coupling effects are expected to have a greater influence over the second 


mode frequencies than they do for the first mode frequencies. 


Figures 23 and 24 give the first and second mode shapes plotted versus the 
NASA measured mode shapes. The data points represent each segment’s deflection for 
the Holzer results and the solid line connects the data for the NASA results. It is 
important to note that the span locations for the models used do not match between the 


Holzer and NASA results. 
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Natural Frequency = 289.5 rad/sec 
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Figure 23. Full-Scale UH-60A Free-Free First Mode Shape. [After: Ref. 15] 
From mid-span to tip the Holzer mode shape compares well with the 
NASA mode shape, but the match is poor from root to mid-span where the mass and 


stiffness is greater. 
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Figure 24. Full-Scale UH-60A Free-Free Second Mode Shape. [After: Ref. 15] 
For an uncoupled solution, the Holzer method compares reasonably with 
the NASA mode shape. These differences are consistent with the errors presented for the 


second natural frequency. 


The Holzer program is valid for the first natural frequency and mode 
shape calculations for a realistic, Free-Free, full-scale, non-uniform rotor blade. The 
Holzer program results for the second mode shape and frequency calculations are an 
acceptable match considering the solution is uncoupled. Use of the Holzer program for 
higher than the first mode of a free-free, non-uniform rotor blade should be confirmed 
with other methods to validate the results. A possible source of error is the fact that the 
Holzer program did not account for shake test boundary conditions. The shake test 


boundary conditions include the effects of the shaker, bungee cords and gravity [Ref. 15]. 
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2. United Technologies Research Center (UTRC) Model-Scale UH-60 
Flex-Free Rotor Blade 


a. Description of UTRC model-blade 


In 1984 UTRC manufactured a four-bladed (1:5.73) model-scale UH-60A 
rotor system developed for wind tunnel testing. The blade radius is 56.22 inches, both 
the radial and feathering axis pitch link offsets are three inches and the chord is 3.643 
inches. The blade aspect ratio is 15.3 and the solidity is 0.0825. The blade is 
dynamically and geometrically scaled. Two different airfoils are used throughout the 
blade span. Figure 25 shows the actual size of the model-scale airfoils. Pressure 


transducers; strain gages and temperature sensors were installed in the blades. [Ref. 14] 


SA 


a) SC 1095 sirfoil. 


CO 


b) SC 1095K8 airfoil. 
Figure 25. Airfoil Sections, Actual Size. [From: Ref. 14] 


The rotor blades are constructed from graphite-epoxy spars, fiberglass 
skin, foam trailing edge pockets, and tungsten counterweights [Ref. 14]. Tables 16 
through 20 in Appendix C, present the structural properties based on the UTRC data [Ref. 
17]. For all the structural property tables the value increase from left to right and 
coincide with the span locations in Table 16. The tension-torsion term is presented in 
Table 19 without the square of rotation speed, which is multiplied internally in the Holzer 
program. The tension-torsion term is already multiplied by the square of the radius of 
gyration. The torsional stiffness is determined by measuring deflection resulting from 
point calibration loads. The instrumented rotor blade is thirty-percent heavier than the 
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scale weight of a full-scale blade without instruments. The added weight caused natural 


frequencies to change by as much as ten-percent. [Ref. 14] 


The root restraint was modeled to be flexible like the actual UH-60A 
model rotor . The value of kgoot is 10,300 (in-lb/rad) and Equation 2.14 results in a root 
torsional deformation boundary condition of 0.0904 rad. [Ref. 14] 


b. Non-Rotating Natural Frequencies and Mode Shapes 


This test case is the first examination of the Flex-Free root boundary 
condition. In addition, this is the first test case to investigate the effects of the tension- 
torsion and “tennis racket” terms on natural frequencies and mode shapes. Figure 26 
shows the graphical solution to the Flex-Free Condition. The root flexibility allows for 
deflections up to 0.0904 rad. The Holzer program assumes the maximum allowed 
deflection will be achieved at the natural frequencies. The odd modes have a positive, 


kproot, and the even modes have a negative, kroot. [Ref. 1] 


O root B.C. met 
— full set of root solutions 
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odd mode naturaNfrequencies 
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Figure 26. Graphical Solution to the Non-Rotating, Non-Uniform, Flex-Free Natural 
Frequency. 
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Table 11 presents the natural frequency solutions for the non-rotating, 
Flex-Free Holzer calculations, the Myklestad based Sikorsky calculations, and rap test 


data. Table 11 also shows the calculated natural frequencies for Free-Free solutions. 


Flex-Free} error | Free-Free| error 
a | % | ate | % 
124.8 


Sikorsiy +866 | 69 | 88 8 
Rap et fe gO — rs 


Table 11. | Natural Frequencies for the Non-Uniform, Non-Rotating, Flex-Free Rotor. [After: 
Ref. 14] 


Results of the Holzer program’s Flex-Free first natural frequency agrees 











better with rap test data than do Sikorsky calculations. The Free-Free first frequency 


result for the Holzer method is in slightly closer agreement with rap test values than the 


Flex-Free value. 


Figures 27, 28 and 29 show first, second and third mode shapes 
respectively. Data points represent individual segment deflections with the line 


connecting the data points. 
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Figure 27. First Mode Non-Rotating Model Rotor Natural Frequency. 
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Figure 28. Second Mode Non-Rotating Model Rotor Natural Frequency. 


Natural Frequency = 2074.5 rad/sec 
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Figure 29. Third Mode Non-Rotating Model Rotor Natural Frequency. 
There are no exact mode shapes to compare the Holzer generated mode 


shapes against. It is noted that the mode shapes do follow expected trends for the first 
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three modes for a Fixed-Free case. The Holzer program results for the non-rotating, 
model-scale UH-60A rotor are valid. 


Cc. Rotating Natural Frequencies and Mode Shapes 


The ultimate purpose for the Holzer method described here is to calculate 
mode shapes for a rotating, non-uniform, Flex-Free rotor blade. This section explores 
results of the Holzer method with and without inclusion of the tension-torsion and “tennis 
racket” terms. The case considered has a rotational speed of 1,442 RPM (151 rad/sec or 
24 Hz). Sikorsky has calculated the natural frequencies for the UTRC blade. Sikorsky 
calculations are based on a coupled Myklestad method. Table 12 shows Holzer method 
and Sikorsky calculated rotating natural frequencies and deviation. There are no 


measured or calculated mode shapes available for comparison with the Holzer generated 


@, dev @) 
Rotating, Q=1442 RPM (Hz) % (Hz) (Hz) 


3.0_| 250.1 | 330.9 
Sikorsky (coupled Myklestad) | 89.1 | | 


Table 12. Natural Frequencies for Non-Uniform, Rotating, Flex-Free Rotor. [After: Ref. 14] 


mode shapes. 














Observe that the first uncoupled Holzer natural frequency is three percent 
less than the Sikorsky coupled value. No comparison solutions are available for the 
second and third modes. Figure 30 shows the differences between the rotating and non- 


rotating mode. The lines connect the deviation data points for each mode shape. 
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Figure 30. Deviation (x1 0° ) Between Rotating and Non-Rotating Mode Shapes. 


The variation between rotating and non-rotating mode shapes is 
insignificant with the largest deviation being less than —0.002 for the third mode (tip 
deflection equals 1.0) as shown in Figure 30. Another way to analyze the effect of 
rotational effects on natural frequency calculations is with a Southwell plot. Figure 31 is 
a Southwell plot displaying the first four elastic mode natural frequencies as they vary 
with rotor speed with Sikorsky coupled data shown in solid lines [Ref. 14]. Circles at 
zero RPM mark the rap test results [Ref. 14]. Multiples of rotor speeds are represented 
by dotted lines that are labeled along the right side of the plot as 1P (one-per revolution) 
through 6P. The Holzer calculated, uncoupled torsion natural frequencies are marked by 


an “x” every five RPM and are connected by a dotted line. 


47 


© Rap Test Data 
—<— Holzer Uncoupled 1st Torsion Mode 
— Sikorsky Coupled Modes 


3rd Flap/Lag 


2nd FlapiLae 


Frequency, (Hz) 





in) 400 600 800 1000 1200 1400 1600 
Rotor Speed, (RPM) 


Figure 31. Southwell Plot of Model UH-60 Blade Modes. [After: Ref. 14] 


The Sikorsky calculations show strong coupling between the first torsion 
and second flap modes. Where this coupling is strongest, between 600 and 1,000 RPM, 
the frequencies of these coupled modes are widely separated from the uncoupled torsion 
mode calculated with the Holzer method. Above 1,000 RPM, however there is good 


agreement between the Sikorsky calculations and the Holzer method. 


Based on the results presented the Holzer method is valid for the Flex- 
Free, rotating, non-uniform rotor blade case. The Holzer method has been demonstrated 
to be of sufficient accuracy to apply to the calculation of mode shapes required for Strain 


Pattern Analysis (SPA) 
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IV. DNW RUN 13.20 OF THE MODEL UH-60 ROTOR BLADE 


A. INTRODUCTION 


To test the Strain Pattern Analysis (SPA) program listed in Appendix B, an 
available database of strain measurements and the rotor blade structural properties must 
be available. The structural properties are used in the Holzer program (Appendix A) to 
generate mode shapes needed for the SPA method. The SPA method can be tested by 
applying the calculated mode shapes and measured torsion moment data to Equations 1.7 
through 1.9. The model UH-60A blade used in the final Holzer test case in Chapter III 
has an extensive database of wind tunnel test results available. The database includes 
torsional strain measurements at six span locations that are sufficient to use as the test 


case for the developed SPA program. [Ref. 14] 


The 1989 test of the UTRC manufactured (1:5.73) scaled UH-60A rotor blade at 
DNW represented the most extensive model rotor test ever conducted at that time. The 
test program is documented in a seven-volume set of reports plus digital data tapes. The 
UTRC model blade is described in Chapter III.C.2. The rotor system is a four-bladed 
(1:5.73) model-scale UH-60A rotor system developed for wind tunnel testing. The blade 
radius is 56.22 inches, both the radial and feathering axis pitch link offsets are three 
inches and the chord is 3.643 inches. The blade aspect ratio is 15.3 and the solidity is 
0.0825. The test measured blade pressures, acoustics, dynamics, and performance during 
numerous test runs including 42 separate level flight conditions, 43 descent conditions 
and 19 hover conditions. Test run number 13.20 was chosen from the level flight 
condition runs based on availability of the data from previous correlation studies. 
Instrumentation installed throughout the four rotor blades included 176 pressure 
transducers, 16 strain gages, eight temperature detectors, and 12 hot film anemometers. 
Each strain gage measured 1,024 samples per rotor revolution. For the purpose of this 
thesis only the torsional strain gage measurements are used. The rotor system was 
mounted on a fully articulated hub and the hub was mounted onto the Army 
Aeroflightdynamics Directorate (AFDD) rotary wing test stand. The test program was a 
joint effort between the AFDD, NASA, UTRC and Sikorsky Aircraft. [Ref. 14] 
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B. DNW FACILITY AND TEST AND RUN 13.20 CONDITIONS 


The German-Dutch wind tunnel, (DNW), is a joint aerospace research laboratory 
run by Germany’s DLR and the Netherlands’ NLR research organizations. The wind 
tunnel complex opened in 1980 and is Europe’s largest and most extensive wind tunnel. 
Test data for the model UH-60A blade was taken in the open-jet configured tunnel with 
maximum flow velocities of up to 262.5 ft/sec (155 knots). The test area dimensions are 


(170 x 98 x 65 ft). Table 13 provides the details for DNW run 13.20. [Ref. 14] 


UH-60A MODEL ROTOR, RUN 13, POINT 20 date: 3/22/89 
213.26 ft/sec 


24 Hz 
1442 RPM 
151 rad/sec 


| 0 | 90 | 180 
0.6324 | 0.8224 


Table 13. .DNW Run 13.20 Data. [After: Ref. 14] 
The advance ratio is 0.301, the advancing blade Mach number is 0.8224 and the average 
Reynolds number is 1,278,729. 
C. DNW MEASURED TORSIONAL MOMENTS 





Three (four-active arm) strain gage bridges are located at five span locations (0.2, 
0.35, 0.5, 0.675 and 0.745). At each of the five locations the first gage measures flapwise 
bending, the second gage measures chordwise bending and the third gage measures 
torsion moments. The sixteenth strain gage is located at 0.875 r/R and measures torsion. 
The torsional moments were determined by applying a calibration matrix from static load 
tests to the strain gage output voltages. Figure 32 presents the torsion moment time 


histories at all six span locations. [Ref. 15] 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 32. Torsion Moment Time History for DNW Run 13.20. 

The torsion moments in Figure 32 have their r/R span locations labeled on their 
curves from 0.2 to 0.875. Each span location’s torsion moment is transformed into 
harmonic terms using a Fourier transform and the first ten terms are used in the SPA 
program. 


D. DNW STATIC DEFLECTION CALCULATION 


Dr. Lorber of UTRC [Ref. 14] calculated the torsional deflection as described in 
the description of the static method in Chapter I. The DNW results were replicated by 
applying the static method to the DNW measured torsion moments and given stiffness 
property of the model-scale blade. Figures 33 to 42 present DNW’s calculated torsional 
deflection compared with the replicated static method calculations. The torsional 
stiffness from Table 17, Appendix C and the measured torsion moments for run 13.20 
were interpolated to a step size of 0.5 r/R. The tip torsion moment was assumed to be 


zero and the torsion moments between the last strain gage at 0.875 1/R and 1.0 were 


al 


assumed to linearly decrease to zero at the tip. The root torsion moment was extrapolated 
to the 0.534 radial axis pitch offset location. [Ref. 14] 
UH-60A Model Rotor, DNW Run 13.20, 4 =0.301 
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Figure 33. Static Method versus DNW Calculated Torsional Deflection at 0.225 r/R. 
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Figure 34. Static Method versus DNW Calculated Torsional Deflection at 0.4 r/R. 
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Figure 35. Static Method versus DNW Calculated Torsional Deflection at 0.55 r/R. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 36. Static Method versus DNW Calculated Torsional Deflection at 0.675 r/R. 
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Figure 37. Static Method versus DNW Calculated Torsional Deflection at 0.775 r/R. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 


Torsional Deflection (deg) 





“O 50 100 150 200 250 300 350 
Azimuth (deg) 


Figure 38. Static Method versus DNW Calculated Torsional Deflection at 0.865 r/R. 
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Figure 39. Static Method versus DNW Calculated Torsional Deflection at 0.92 r/R. 
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Figure 40. Static Method versus DNW Calculated Torsional Deflection at 0.945 r/R. 


Torsional Deflection (deg) 





“o 50 100 150 200 250 300 350 
Azimuth (deg) 


Figure 41. Static Method versus DNW Calculated Torsional Deflection at 0.965 r/R. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 


Torsional Deflection (deg) 





“o 50 100 150 200 250 300 350 
Azimuth (deg) 


Figure 42. Static Method versus DNW Calculated Torsional Deflection at 0.99 r/R. 

The replicated torsional deflections from the 0.55 span location to 0.99 r/R are 
excellent. Minor differences in Figures 35 through 42 could be the result of differences 
in interpolation and rounding. The differences are more apparent for the span locations 
from 0.225 to 0.4. The differences near the root are due to different integration step sizes 
between the replicated static method and the DNW calculations. The DNW static 


deflection results will be used to evaluate SPA method deflections. 


af 
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V. APPLICATION OF STRAIN PATTERN ANALYSIS 


A. INTRODUCTION 


A comparison between the developed SPA method and the DNW static method 
will be conducted for the torsional deformation results of the UTRC model UH-60A rotor 
blade during run 13.20 of the DNW test presented in Chapter IV. The SPA method is 
applied to the Holzer program generated mode shapes as presented in Chapter II for the 
UTRC blade. The application of the SPA method is done with the SPA program listed in 
Appendix B. The SPA method follows steps one through three as presented in Chapter I 
with calculations based on Equations 1.7 through 1.9. 


B. SPA STEP ONE 


The SPA vector, {q}, is derived for the first 21 harmonics using Equation 1.7 


restated below. 


{q} = (TPT ET] {M} (1.7) 
Each of the first 21 measured harmonic terms of the measured torsion moment is 
represented as a vector {M}. The calculated torsion moment mode shape matrix, [T], is 
an output of the Holzer program. There were six strain gage locations that were used to 
measure torsion moments for the UTRC blade. From the SPA theory in Chapter I, the 
maximum number of modes for use in the SPA method (columns of the [T] matrix) is 
five, one less than strain gage locations. The SPA program will be run five separate 
times to study the effect the number of modes used in the SPA method has on the results. 
The five separate SPA solutions are for each number of mode shapes used from one 


mode shape to the maximum of five mode shapes. 


In Appendix D, Figures 44 through 53 show the solution for the SPA vector, {q}, 
for all 21 harmonic terms for the five-mode SPA solution. The lines connect the data 
points for each harmonic. The value of the SPA vector term for each mode determines 
that mode’s contribution to the least squares fit of the measured torsion moment and 
ultimately the deflection solution. The contribution of each mode decreases for higher 


modes. The SPA vectors converge to zero as they approach the fifth mode. 


59 


C. SPA STEP TWO 


To confirm the accuracy of the SPA vector {q} Equation 1.8, given below, is 


used. 


{M} =[T]{q} (1.8) 
The SPA vector is multiplied with the calculated torsion moment matrices, [T], from the 
five runs of the Holzer program. The result of Equation 1.8 is the SPA method’s least 
square fit of the measured torsion moment data for the first through fifth mode SPA 
vector. The estimated torsion moment is plotted against the measured torsion moment on 
Figures 54 through 68 in Appendix D for each of the first 21 harmonics for the first 
through fifth mode solution. A “+” marks the measured torsion moment values and the 
dotted line connects the measured moment data points. A diamond marks the SPA 
calculated torsion moment values resulting from Equation 1.8, and the solid line connects 
the SPA data points. We observe that the results are poor for the one-mode SPA torsion 
moment estimation on Figures 54 through 56. The two-mode SPA torsion moments in 
Figures 57 through 59 are a great improvement over the one-mode solutions, but are still 
a poor fit to the measured data. In Figures 60 through 62 the three-mode SPA torsion 
moments are shown to be the first SPA solutions to accurately fit the measured data. The 
three-mode SPA results significantly miss the peaks and valleys of the measured data. 
The four-mode SPA results shown in Figures 63 through 65 are an improvement upon the 
three-mode solution, but still significantly miss some peaks and valleys of the measured 
data. Figures 66 through 68 show the five-mode SPA results are clearly the best fit of the 
measured data. The SPA method will continue to be applied using only the five-mode 


SPA vectors. 


Figures 69 through 74 show the time histories of the five-mode SPA torsion 
moment compared to the measured torsion moment for each strain gage location. The 
measured moment dotted line clearly shows where the discrepancies between the SPA 
results and measured values are located. The five-mode SPA fits the measured torsion 
moment time histories well, but the fit is not perfect. Higher harmonic influences may be 


the cause of the mismatch between the SPA moments and the measured moment values. 
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D. SPA STEP THREE 


The final step of the SPA method applies Equation 1.9, shown below, to estimate 
the deformation, {@}, using the torsional deformation matrix, [D], multiplied by the five- 


mode SPA vector, {q} derived in step one and validated in step two. 


{o} =[Dlig} (1.9) 


The multiplication is done twenty-one times, once per harmonic. Figures 75 to 84 plot 
the time histories of the SPA torsional deflection against the DNW static method results 
from Chapter IV for each reported span location. It is important to note that there is no 
exact solution available for this test case. The deviations between the DNW results and 
the SPA results are not an indication of either method’s inaccuracy. The overall 
similarities between the two methods establish confidence in both methods’ validity. 
Optical method results are unavailable for comparison. Both methods calculate a similar 
deflection pattern at each span location. The maximum deviation between the two 
methods is approximately 1.1 degrees and occurs at approximately 172 degrees azimuth. 
There are also relatively large deviations in the region between 60 and 90 degrees 
azimuth. The differences between the two methods in the noted locations are large as a 
percentage of total deflection. Tables 14 presents the maximum, minimum and mean 


deflection values for each span location. 


UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 


FR [max] min [mean [max] min_| 
: 2 


m 
0.225 -0.33 
0.400 | 0.83 | 347 | -0.99 | 0.74 | -240 | 
0.55 -0.87 
0.67 -0.95 


0.86 
0.92 0.88 
0.94 0.85 
0.96 0.83 
0990 [2.61 | 554 | -099 | 249 | -466 | -0.82 


Table 14. © Maximum and Minimum Torsional Deflection by Span Location. 
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0.77 





The mismatch between the static method and SPA method may be a result of 
higher harmonics influencing the static method result versus the SPA method that only 
includes harmonics up to the 10th sine term. An attempt was made to include harmonics 
through the 20th sine and no significant change in the SPA result was identified. Figure 
43 shows the mean values of SPA and DNW deflection from Table 14 plotted versus 
span location. The lines connect the data points for the SPA deflections and the DNW 
deflections. It is important to only use this figure in the context of the actual deflection 
time histories or the significance of torsional deformations around the entire rotor 
azimuth will be overlooked. Comparing the maximum and minimum values in Table 14 
with Figure 43 is a good indication that the average deflection values do not give a 


complete picture of torsional deformation. 
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Figure 43. Five-Mode SPA versus DNW Average Deflection by span location. 
Figure 43 shows how the overall pattern of the SPA deflection is similar to the 
DNW deflections, but the SPA deflections are offset by approximately three-tenths of a 
degree from the root to 0.865 r/R. The average (steady) values for torsional deflection 


are sensitive to the number of modes used in the solution. To analyze the differences 
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using the first ten harmonics of the SPA estimated deflections and the DNW static 
method deflection calculations Figures 85 to 94 present linear regression plots for each 
reported span location. Figures 85 to 94 plot the SPA deflection on the y-axis against the 
DNW deflection on the x-axis. Any deviation from the diagonal axis signifies a 
deviation between the two techniques. Figures 85 to 94 show that many of the harmonic 
terms are huddled near the zero axes and are therefore insignificant contributors to the 
deflection solution. The significant terms are marked with the harmonic term they 
represent. Close to the root only the steady, first cosine, first sine and third sine terms 
stand out as strong influences to the deflection solution. Throughout the entire span these 
four harmonic terms continue to dominate the behavior of torsional deformation. As the 
tip is approached six additional harmonic terms significantly influence the torsional 
deflection including the fourth cosine, third cosine, seventh sine, second sine, fifth cosine 


and fourth sine terms. 


Along the bottom of the regression plots the y-intercept (b), correlation 
coefficient, and slope are displayed for the line that best fits the data (the line is not 
plotted). The correlation coefficient is the statistical measure of the degree of association 
between the x-axis (DNW) and y-axis (SPA) data [Ref. 8]. The correlation coefficient 
equals the covariance between the DNW and SPA data divided by the product of their 
standard deviations [Ref. 8]. For a perfect match between SPA and DNW deflection the 
y-intercept would be equal to zero, the correlation coefficient would be equal to one, and 
the slope would be equal to one. Figures 95 through 97 plot the values of the y-intercept, 
correlation coefficient and slope versus span. In Figure 95, the y-intercept remains 
relatively constant just below the ideal zero. The correlation coefficient converges 
towards one as the span locations approach the tip. This indicates a better match between 
the two methods from 0.865 to 0.99 span. The good correlation near the tip can also be 
seen in the time history plots discussed previously. The slope near the root is well off the 
ideal value of one. Near the root the deflections are small due to high stiffness values and 
the harmonic terms cluster near zero and lose statistical significance. The harmonic 
terms increase in significance as the stiffness decreases towards the tip and the regression 


plots and the slope values become a significant indication of the better agreement 
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between the SPA and DNW results toward the tip. As the span locations increase 


towards the tip the slope values converge towards the ideal value of one. 


The similar torsional deformation results for the five-mode SPA method and the 
static method verify the applicability of each method in estimating torsional deflection 
from measured strain data. No further analysis can be accomplished without optical 
measurements to compare with the SPA results. 

E. EVALUATE THE ROBUSTNESS OF THE SPA METHOD 

1. Effect of Changes in Mass and Stiffness Properties 


The SPA method is dependent upon the measured structural properties of stiffness 
and mass. Errors in the measurement of the structural properties manifest themselves in 
inaccurate mode shapes and ultimately incorrect deflection estimates. The Holzer 
program includes an option to randomly vary the mass and stiffness properties plus or 
minus a chosen percentage at each span location to test the SPA method’s robustness to 
structural property inaccuracies. The five-mode SPA test case for the UTRC model UH- 
60A blade was recalculated using a random five percent variance in structural properties. 
Figures 98 and 99 present the deflection results using varied structural properties on time 
history plots for the 0.4 and 0.775 radial stations compared to the original five-mode 
deflection result. The result of varying mass and stiffness properties appears to be 
insignificant to the validity of the SPA deflection result. 


2. Effect of Number of Measurement Stations 


A strain gage failure during a test program can also affect the ability to estimate 
deflections using the SPA method. The SPA program allows for the elimination of all 
data from a failed strain gage. To test the robustness of the SPA method to eliminating 
one strain gage’s measurements from the torsion moment data another run of the Holzer 
and SPA programs was conducted where the data from the strain gage at 0.642 was 


eliminated. 


Figures 100 and 101 show the time histories for the 0.4 and 0.775 radial stations 
compared with the original five-mode SPA result. The maximum number of modes 
allowed to fit measured data from five strain gages is four modes. With only six strain 


gages in the DNW test of the UH-60A model blade the loss of any one strain gage has a 
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significant impact on the final SPA result. The loss of one strain gage has a two-fold 
impact on the SPA technique. First, there is less available data to fit the sets of calculated 
modes to. Second, the SPA solution process must reduce the maximum number of modes 
used to estimate the deflection. 


3. Effect of Measurement Error 


Short of a strain gage failure it is normal to have errors in experimental 
measurements of rotor blades [Ref. 3]. Measurement errors include calibration error, 
static error, and random error [Ref. 3]. The SPA program was run for five modes using a 
random error from minus five to plus five percent at each of the six strain gage locations 
on the DNW blade. Figures 102 and 103 show the time history plots for the deflection 
results at 0.4 and 0.775 radial stations compared with the original five-mode SPA result. 
The solution for each harmonic term changed negligibly from the original five-mode SPA 
test case. The method is insensitive to the randomly occurring measurement error of five 


percent. 
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VI. SUMMARY 


A. CONCLUSIONS 


The Strain Pattern Analysis (SPA) method has been developed and tested for 
torsional deformation using data from the model-scale UH-60A rotor blade wind tunnel 
test conducted at DNW. The Holzer method was used to provide the SPA method with 
the calculated torsion moment mode shapes and torsional deflection mode shapes for the 
rotor blade. The measured torsion moments are from wind tunnel tests conducted by the 
Army, NASA, UTRC, and Sikorsky at DNW with a (1:5.73) model-scale UH-60A rotor 
blade with an advance ratio of 0.301, an advancing tip Mach number of 0.8224 and an 
average Reynolds number of 1,278,729. The SPA method estimates torsional deflections 
in two steps. First, the mode shapes of a rotating blade in a vacuum are calculated. 
Second, the modal amplitudes are identified using a least squares fit to the measured data. 
The mode shape calculation only has to be done once unless more than one rotor speed is 


evaluated. 


The SPA method torsional deflections compare well with the overall trend and 
range of static method integrated deflections from the DNW database. The SPA method 
predicts slightly larger deflections than the static method. Differences between the two 
methods of estimating torsional deflection are expected. No additional validation of the 
SPA results can be done without results from optical methods of measuring torsional 


deformation. 


Possible sources of error in the application of the SPA method include the Holzer 
method’s approximate solution to the equation of motion, inaccurate strain gage 
measurements and inaccurate rotor blade structural properties. The Holzer calculated 
mode shapes are based on the developed uncoupled rotor blade equation of motion for 
torsion. Coupling effects could change the mode shapes and the final SPA deflection 
estimate. The possibility of a faulty strain gage was tested and the SPA program proved 
to be sensitive to the decrease in available strain gages for the six strain-gage UH-60A 
blade. The SPA method was demonstrated to be tolerant to variations in blade properties 


and measurement errors. 
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B. RECOMMENDATIONS FOR FURTHER RESEARCH 


Further areas of study should include the evaluation of the SPA method for other 
test runs from the DNW database for the model UH-60A rotor blade. A complete 
evaluation of the SPA program developed in this thesis would require data from test 
programs with available optical measurement results to compare to. Improvements can 
also be made in the calculation of rotor blade mode shapes. The Holzer program 
developed in this thesis can be improved by using the coupled equations of motion and 
developing the Myklestad method for determining natural frequencies and mode shapes 
as presented in [Ref. 9] and [Ref. 18], in which the fully coupled Myklestad method is 
developed for specific application to rotary wing aircraft. The sensitivity towards a 
decrease in modes caused by a simulated strain gage failure should be investigated 
further by using test cases with more measurement stations. The higher mode’s 
contribution to the SPA deflection estimate is expected to become increasingly 
insignificant and therefore the method should prove tolerant to the loss of a strain gage 
for tests involving a greater number of measurement stations. Further research could 
determine the minimum number of measurement stations required for the SPA method to 
be insensitive to the loss of a strain gage and a decrease in the number of contributing 


modes. 
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APPENDIX A. HOLZER PROGRAM 


A. INSTRUCTIONS FOR RUNNING THE HOLZER PROGRAM 
1. Holzer Program Set Up 


To start the Holzer program type “SPAholzer” in the MATLAB™ command 
window. The program queries the user for the information necessary to calculate the 
torsional mode shapes. After you select to enter your own test case or to use one of the 
built-in test cases (line 8) you will be asked for the number of modes to solve for (line 
20). The number you select will determine how many mode shape plots are displayed 
and how many modes are used in the SPA solution. Choose zero modes if you choose to 
solve for natural frequencies over a range of rotational speeds or to individually choose 
modes for plotting after the natural frequencies are found. For brevity the program listing 
presented does not list the structural properties for the built-in test cases. The next option 


selects the Fixed or Flex-Free versus Free-Free boundary condition (line 30). 


If you opt to enter your own data you will first be asked for the number of stations 
you have structural properties for (line 32). If you choose to enter your own data it is 
important to be consistent with units. The program asks for (lb-ft-sec”) for the torsional 
mass moment of inertia (line 41), (fi”-Ib) for torsional stiffness (line 43) and (ft) for the 
distance between segments (line 46). In all cases the first segment is the one closest to 
the root. Table 15 gives the input units you may use and the corresponding output 
variable units. The input and output units for rotational speed and frequency will always 


be rad/sec. 


[1 Gi] _ | Torsion Moment [Deflection 


Table 15. Input and Output Units of Measure. 





The built-in plot labeling will always give the English units in feet and is easily 
changed using the option to change axes properties from the edit option on the figure 
toolbar. If you will be using your own test case you must have matrices of structural data 


ready to paste into the command window. When the program queries for a specific 
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variable copy the corresponding matrix from its location, then on the command window 
type a matrix symbol “[ ]” and paste your matrix in the middle of the brackets. Each 
structural property matrix must have the same number of terms, but the matrices can have 
different numbers of rows and columns. If you will be using centrifugal terms you must 


prepare the blade angle prior to deformation, /, in degrees for each station (line 105) 


and the tension-torsion term (lines 33-38) in the form, ((1/2)(mass/Ar)(R? —r?) *k,’). 
In this form the Holzer program will multiply by the appropriate rotational speed to 


include the stiffening effect of the tension-torsion term. 


Whether you enter your own data or use a built-in test case you will be presented 
with a number of options prior to solving for natural frequencies and mode shapes. You 
may opt to not include the “tennis racket” term in your calculations (line 103). If your 
structural properties do not extend all the way to the tip (1.0) location you can apply the 
properties from the closest station to the tip to create a final station where deflection will 
equal one (line 40 and lines 195-207). This option forces the boundary condition of one 
radian tip deflection to be applied at the tip instead of the last location available. These 
options are designed to allow the effects of individual terms to be evaluated on the final 


result. 


The program will ask for the measurement stations for test data you have for use 
with the SPA program (lines 49 & 51). This information is used to interpolate the 
calculated torsion moment mode shape to the same span locations as the measured data to 
allow for the SPA vector to be calculated (lines 370-378) for Fixed & Flex-Free, (lines 
430-439) for Free-Free. Enter the locations between “[ ]” symbols in non-dimensional 
r/R form. After choosing your rotational speed or range of rotational speeds (lines 126- 
137) you will be asked for test frequency parameters (lines 139-150). The most useful 
method of solving natural frequencies is to enter a range of “test” frequencies from zero 
(lower limit) to the an applicable upper limit that will allow for your desired number of 
modes to be found. I suggest using the highest mode you wish to solve for times 1,000 as 
an upper limit (i.e.: if you want five modes, enter 5,000 as your upper limit). The Holzer 
program will locate all the natural frequencies between your lower limit and upper limit 


by testing each frequency from zero to the upper limit at the entered step size. Your last 
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options are to include a random variation to the mass and stiffness properties to allow for 


a test of program robustness (lines 156-181). 


The first step of the Holzer method is Equation 2.10 on line 211. Step two is on 
line 212 using Equation 2.21. Step Three is on line 218 using Equation 2.23. Step Four 
is on line 267 with Equation 2.21. The program loops as many times as required to 
repeat steps three and four until the root solution for each tested frequency is calculated. 
When each test frequency root solution is completed the program will execute step six to 
find the natural frequencies. Step six uses the bisection method to find where the root 
deflection equals zero for the Fixed-Free condition or the allowed deflection for the Flex- 
Free condition on lines 221 to 263. For the Free-Free condition the bisection method 
searches for where the root torsion moment equals zero on lines 270 to 292. 


2. Holzer Program Operation 


As the Holzer program finds natural frequencies it will display the frequency in 
rad/sec, the mode number, and the rotational frequency in rad/sec on the screen. The 
program will display a graphical solution plot (lines 306-320/322-329) to the natural 
frequency when it completes the full range you entered. The graphical solution can be 
used to confirm no natural frequencies were missed. It is possible to miss a higher mode 
when the step size is larger than one. If a mode is missed you must start over and use a 
smaller step size. You are given the option to delete a measurement station location to 
test the effects of a bad strain gage (lines 350-361/410-429). The program will now 
display the deflection mode shapes for only the number you selected (lines 384-401/447- 
459). If you chose zero for number of modes you will have the option to individually 
enter frequencies and rotational speed to plot mode shapes for (line 464). The outputs 
from the Holzer program are listed in the next section. 


B. HOLZER PROGRAM VARIABLES 


DD (n x m) Torsional Deformation Mode Shapes (rad) where each row 
corresponds to the stations structural properties were provided for and each column is a 


separate mode (line 346 for Fixed/Flex-Free) and (line 406 for Free-Free). 


yal 


GG (n x m) Torsion Moment Mode Shapes ( ) where each row corresponds to 
the stations where structural properties were provided and each column is a separate 


mode (line 347 for Fixed/Flex-Free) and (line 407 for Free-Free). 


G (g x m) Torsion Moment Mode Shapes ( ) interpolated to measurement 
station locations for rows and each column is a separate mode (line 377 Fixed/Flex) and 


(line 438 Free-Free). 


WNF# WNF1 is the first natural frequency (rad/sec) in order of increasing 
rotational speed (rad/sec) and WNF2 is for the second natural frequency (line 240 & 
245). 


OMEGA# The rotational speeds (rad/sec) that correspond to WNF1,WNF2 (line 241 
&246). 


rootfreq The natural frequencies (rad/sec) in mode order if a single rotational speed 


was selected (line 234 Fixed/Flex) and (line 282 Free-Free). 


x(n) The non-dimensional calculation stations where structural properties were 


provide, increases by one station if you forced the tip boundary condition (line 98). 


span(g) The non-dimensional measurement stations (line 51). 
C HOLZER PROGRAM LISTING 

SHolzer's method for ROTATING Torsion Mode approximation 
2 le 


1 

end 
3 clear 

4 eliminate=2; 
5 

6 

4 





screate the 1x1024 matrix of azimuth locations for DNW cases 
primarily 
NEWAZIMUTH=[0: (360/57.3) /1024: ((360/57.3) -(360/57.3)/1024) ]; 
CNEWAZ IMUTH=NEWAZIMUTH (:) ; 
8 choice=input ('enter 1 to input your own blade properties \n or 2 
for 126 segment DNW\n or 3 for Full-Scale UH-60A data\n or 4 for the 
Scanlan Wing example\n or 5 for the 49 segment uniform, fixed test 
case\n'); 














10 if choice== 

11 col=1; 

12 end 

13 yl=1;%These variables are needed for WNF1, OMEGA1 for southwell 
plots 

14 y2=1; 

15 disp('for all inputs you must use a full matrix') 

16 disp('add zeros if required even if it exceeds your number of 
nodes') 

dey disp('the total number of elements in each matrix must be the same 
size') 


18 disp('but rows & columns do not have to be the same length") 
19 disp(' '); 





ip 


20 modenumber=input ('enter 0 to individually select (also for range of 
omegas) \n or enter desired # to preset how many natural frequencies to 
do mode shape plot and/or SPA analysis for\n'); 

22 modec=modenumber; 

23 CTLSTIFF=input ('For Flex-Fr nter the root torsion displacement\n 
(zero for Fixed-Free and Free-Free, .0904 rad for DNW)\n'); 

24 if CTLSTIFF~=0 

25 mult=2; 
26 end 
27 if CTLSTIFF== 
28 mult=1; 









































29 end 

30 choice4=input (‘enter 1 if this is a Flex-Free or Fixed-Free 
boundary condition or zero for Free-Free\n'); 

31 if choice== 

32 nodes=input ('enter the number of segments\n to evaluate'); 

33 tension=input ('enter 1 to enter your own tension-torsion term\n 
(without rotation speed) or zero to not include\n'); 

34 if tension~=1 

35 TENS=zeros (nodes,1); 

36 end 

37 if tension== 

38 TENS=input ('enter the vector of tension-torsion terms\n') 

39 end 

40 force=input ('enter 1 to force the tip boundary condition or zero to 





apply deflection of one radian at your last station'); 

41 CITHETA=input (‘enter your Torsional Mass Moment of Inertia (lb-ft- 
s*2) matrix, \n'); 
42 CITHETArow CITHETAcol]=size(CITHETA) ; 
43 CGJ=input (‘enter your NON-ZERO GJ (ft*2-l1b) matrix,\n ensure same 


number of elements as segments chosen (theoretical fixed segment gets 
arbitrary non-zero value) \n'); 
























































45 CGJrow CGJcol]=size(CGJ); 

46 Cdeltax=input ('enter your delta r (feet) matrix,\n ensure same 
number of elements as segments chosen, (theoretical fixed segment=0)\n"'); 
48 Cdeltaxrow Cdeltaxcol]=size(Cdeltax) ; 

49 experiment=input ('Enter 1 if you have experimental measured data to 
use with SPA\n or zero if you only want modes calculated\n'); 

50 if experiment==1 

51 span=input ('enter the r/R Strain Gage locations=\n'); 

52 end 

33 if experiment~=1 

54 span=[0:0.1:1]; 

55 end 

56 The following loops are to make column vectors out of all matrices 
entered$$%%%S%S%SSSSS%S 

57 i=l; 

58 j=1; 

59 k=1; 

60 for i=l:nodes 

61 if jJ==CITHETAco1+1 

62 j=1; 

63 k=k+1; 

64 end 

65 ITHETA (i, 1)=CITHETA(k, 4); 

66 i=itl1; 

67 9 j=jtl; 

68 end 

69 i=1; 

70 j=1; 

71 k=1; 

72 for i=l:nodes 

73 if j==CGJcol+1 

74 j=1; 


fe: 


75. k=k+1; 
76 end 
77 GJ(i,1)=CGI(k,4); 
78 i=itl1; 
79 j=5t1; 














80 end 

81 i=1; 

82 j=1; 

83 k=1; 

84 for i=1l:nodes 

85 if j==Cdeltaxcoltl 

86 j=1; 

87 k=k+1; 

88 end 

89 deltax(i,1)=Cdeltax(k,j); 

90 Lets 

91 j=jtl; 

92 end 

93 sumdeltax(1)=0; 

94 for count=2:1:nodes 

95 sumdeltax (count) =deltax (count) +sumdeltax(count-—1); 
96 count=countt1; 

97 end 

98 x=sumdeltax/sumdeltax (nodes) ; 
99 end 

100 SNOW COS2Beta For all cases except DNW which has it built-in 
101 if choice~=0 


if 

102 if choice~=2 

103 choice5=input (‘enter 0 to enter a Beta (collectivetcyclictbuilt-—in) 
matrix or 1 to not include the tennis-racket term\n'); 

104 if choice5== 

105 THETA=input ('enter the Beta sectional blade angle matrix in 

s\n') 




































































106 [THETArow THETAcol]=size(THETA) ; 
107 i=l; 

108 j=1; 

109 k=1; 

110 for i=l:nodes 

111 aif j==THETAcol1+1 

112 4=1; 

113 k=k+1; 

114 end 

115 CTHETA(i,1)=THETA(k,j4); 

116 i=it+l; 

117) j=5+1; 

118 end 

119 COS2THETA=cos (2*CTHETA/57.3); 
120 end 

121 if choice5== 

122 COS2THETA=[ones (nodes,1)]; 
123 end 

124 end 

125 end 

126 choice3=input ('enter 0 to set range or enter 1 to input specific 





omega (rad/sec) \n'); 

127 if choice3== 

128 omega=input ('enter the blade rotational speed you want evaluated 
(rad/sec)\n'); 

129 loweromega=0; 

130 upperomega=0; 

131 end 

132 if choice3== 

133 loweromega=input ('enter lower limit of omega (rad/sec)\n'); 
134 upperomega=input ('enter upper limit of omega (rad/sec)\n'); 
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135 %omega is rad/sec NOT RPM!!! 

136 omega=[loweromega:5:upperomega]; 

137 end 

138 %w is rad/sec 

139 choice2=input ('enter 0 to set range of frequencies or 1 to input a 
specific frequency (rad/sec)\n'); 

















140 if choice2== 

141 w=input('enter frequency to check in rad/sec\n'); 

142 stepsize=1; 

143 upper=0; 

144 lower=0; 

145 end 

146 if choice2== 

147 lower=input ('enter lower limit of frequency search(ensur nough 
range to support # of modes you want solved for) \n'); 

148 upper=input('enter upper limit of frequency search\n'); 

149 stepsize=input ('enter the stepsize between frequency values\n start 








with five, but if a root is missed decrease the stepsize\n'); 

150 w=[lower:stepsize:upper]; 

151 end 

152 zz=1; Sneed this for next line 

153 alternate=(-1)*(zz+1); %Need this for the first mode pass thru 
bisector 





154 ZAAGAFLFA8ALAAAAAAAAALAAAAAAAAAAAAAAAAAAAAAAAAAAAAD 

0000000000000 0000000000000000000000000000000 700 
155 %S%SDesigned to loop as long as desired for specific natl freq 
testing 


156 SThis is to test robustness 

157 mass=input ('enter the max percentage you wish to randomly (+/- from 
zero to max) change TENSION & ITHETA by\n to test robustness of 
solution to mass changes\n'); 

159 stiffness=input ('enter the max percentage you wish to randomly (+/- 
from zero to max)change GJ by\n to test robustness of solution to 
stiffness changes\n'); 

161 if mass~=0 

162 for n=l1:nodes 

163 randomass=randperm (mass) ; 

164 randomsign=randperm (2); 

165 sign=(-1)*randomsign (1); 

166 permass=1+sign*randomass (1) /100; 

167 TENS (n)=TENS (n) *permass; 

168 ITHETA(n)=ITHETA(n) *permass; 

169 n=ntl1; 

170 end 

171 end 

172 if stiffness~=0 

173 for n=l1:nodes 

174 randomstiff=randperm(stiffness) ; 

175 randomsign=randperm (2); 

176 sign=(-1)*randomsign (1); 

177 perstiff=l1+sign*randomstiff(1)/100; 

178 GJ(n)=GJ(n)*perstiff; 

179 n=nt1; 

180 end 

181 end 

182 ccontinue=1; 

183 mode=0; 

184 while ccontinue== 

185 for j=1: ((upperomega-loweromega) /5+1) 

186 if mode== 























187 disp('Natural Frequency(rad/sec) Mode Rotational Frequency 
(rad/sec) ') 

188 end 

189 mc=1; 

190 zz=1; 


vis: 

































































































































































ETA (nodes) *theta(n,i,j); 











from previous theta-STEP THREE 








only at root for each tested 








TROL STIFFNESS%% 





ole 





sign*CTLSTIFF) * (theta(e, (i-1), 4) 


Ils 


191 alternate=1; 

192 for i=1: ((upper-lower) /stepsizetl) 

193 STHE HOLZER METHOD 

194 n=nodes;%SThis will be overwritten if force tip option used 
195 if force== 

196 if x(nodes) ~=1 

197 n=nodes+1; 

198 x(n)=1.0; 

199 deltax(n)=x(n)-x(n-1); 

200 %Assume strux properties for point nodes+1l is same 
201 %as properties at point nodes 

202 GdJ(n)=GJ(n-1); 

203 ITHETA(n)=ITHETA(n-1); 

204 TENS (n)=TENS (n-1); 

205 COS2THETA (n)=COS2THETA (n-1); 

206 end 

207 end 

208 %For the Torsion to the right of the last segment 
209 T(nt1,i,3)=0;%A Holzer Method Boundary Condition 
210 %For the deflection at the last segment 

211 theta(n,1i,j)=1;S5STEP ONE EQN 2.10 

212 %For the Torsion at the last segment-STEP TWO EQN 2.21 
213) T(n,1i,3)=T(nt+1,i, 3) +((w(i)) *2- 

(omega (j))*2* (COS2THETA (nodes,1)))*ITH 

214 Nodes=n; 

215 p=n-1; 

216 for m=1:p; 

217 e=n-1; 

218 %this subtracts previous delta theta 

EQN 2.23 

219 theta(e,i,j)=theta(n,1i,j)- 

T(n,i, Jj) * (deltax(n))/ (TENS (n) * (omega (j)) *2+GJI(n)); 
220 %SSTEP 6: Find the Natl Freq's is done 
frequency 

221 %BISECTION METHOD WITH OFFSET FOR CON 

222 sign=(-1)*(mc+1); 

223 if choice4==1 %fixed free 

224 if e==1 SRoot location 

225 if choice2==0 Srange of w 

226 if i>=2 %$i=1 is at lower limit freg 

227 SThis simply finds where sign changes 

228 if alternate>0 

229 signchange=(theta(e,i,j) 

Sign*CTLSTIFF) ; 

230 aif signchange<0 

231 twotheta=[theta(e,i,j) theta(e, (i-1) 

232 twofreq=[w(i) w(i-1)]; 

233 coeff=polyfit (twofreq,twotheta,1); 

234 rootfreq (mc) =(sign*CTLSTIFF-coeff (2) )/coeff (1); 
235 yaxis(mc)=sign*CTLSTIFF; 

236 [rootfregq(mc) mc omega (j) ] 

237 SOL(zz,:)=[rootfregq(mc),omega(j)]; 

238 %for southwell plots 

239 if mc== 

240 WNF1(y1)=rootfreq(mc) ; 

241 OMEGAI1 (y1)=omega(j); 

242 yl=yl1+l1; 

243 end 

244 if mc== 

245 WNF2 (y2)=rootfreq(mc) ; 

246 OMEGA2 (y2)=omega(j) ; 

247 y2=y2+1; 

248 end 

249 mc=mctl1; 
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250 


endsfor if signchange<0 











































































































251 end%sfor if alternate>0 

252 if (theta(e,i,j)-sign*CTLSTIFF) * (theta ( sign*CTLSTIFF) <0 
293°  22=224¢ 13 

254 alternate=(-1)%(zz+1);%for an even zz we get an odd alternate 
255 end 

256 if CTLSTIFF==0 %Sdon't want to alternate solutions 

257 alternate=1; 

258 end 

259 end%sfor if i>=2 

260 endsfor if choice2==0 Srange of w 

261 end%sfor if e== 

262 end%sfor if choice4==1 %fixed free 

263 %%$%SSSBISECTION ROOT FINDER ENDS HERE FOR FIX FE %SS%S%% 
264 %This is the total Torque, 

265 %SSTEP 4: EQN 2.21 on line 267 

266 %(STEP 5: Repeats Steps 3 & 4 until at root for each tested 
frequency) 

267 T(e,i,j)=T(n,1, 3)+((w(i)) *2- 

(omega (j)) *2* (COS2THETA (e,1))) *ITHETA(e) *theta(e,i, Jj); 

268 n=n-1; 

269 m=mt1; 

270 %SSSBISECTION METHOD WITHOUT OFFSET%%%S%%% 

271 %FOR FREE FREE, EP 6 for when at root for each tested frequency 
272 if choice4==0 %Sfr fr 

273 if e== 

274 if choice2==0 Srange of w 

275 aif i>=2 

276 This simply finds where sign changes 

277 signchange=(T(e,i,j))* (Te, (1-1), 

278 if signchange<0 

279 twotheta=[T(e,1i, 3) Tle, (1-1), 3)]; 

280 twofreq=[w(i) w(i-1)]; 

281 coeff=polyfit (twofreq,twotheta,1); 

282 rootfreq(mc) =(-coeff (2))/coeff (1); 

283 yaxis (mc) =0; 

284 [rootfregq(mc) mc omega (j) ] 
285 SOL(zz,:)=[rootfregq(mc),omega(j)]; 

286 mc=mctl; 

287 end 

288 end 

289 end 

290 end 

291 end 

292 %SSBISECTION ROOT FINDER ENDS HERE FOR FRE 5555S S5S% 
293 This end is for the m=1:p loop 

294 end 

295 i=itl;SNEXT FREQ 

296 %This end is for the range of frequencies loop 

297 end 

298 %This loop works with numerous frequencies and one omega 
299 if choice2== 

300 if choice3== 

301 xl=w; 

302 yl=zeros(size(x 

303 yll(1,1:size(xl -CTLSTIFF; 

304 yul(1,1:size(xl +CTLSTIFF; 

305 if choice4== 

306 figure 

307 plot (rootfreq, yaxis,'s',w,theta(1,:,1),'k-') 

308 line(xl,yl) 

309 line(xl,yll) 

310 line(xl,yul) 

311 if CTLSTIFF~=0 




























































































di 















































312 text(10, (+CTLSTIFF+0.1),'odd mode natural frequencies') 
313 text (10, (-CTLSTIFF-0.1),'even mode natural frequencies') 
314 end 

315 axis([0 upper -1.5 1.5]) 

316 xlabel('Frequency (rad/sec) ') 

317 ylabel(['Torsinal Deformation @ Root (rad)',' ,r/R = 
",num2str(x(1))]) 

318 title('Graphical Solution to Natural Frequency (Fixed-Free) ') 
319 legend('root B.C. met','full set of root solutions'); 
320 end 

321 if choice4== 

322 figure 

323 plot (rootfreq, yaxis,'s',w,T(1,:,1)/1000, 'k-') 





324 line(xl,yl) 

325 xlabel('Frequency (rad/sec) ') 

326 ylabel('Total Torsion Moment (Ft-Lbs) at segment 1 (*10-3)') 
327 title('Graphical Solution to Natural Frequency (Free-Free) ') 
328 legend('root B.C. met','full set of root solutions'); 

329 end 

330 end 

331 end 

332 %if checking a range of omegas this will work 

333 if choice3== 





























334 4 4=j+1;3NEXT OMEGA 
335 end 
336 %This end is for the loop for a range of omega(j) 
337 end 
338 SSSSSESSEEEEEEEEESESESESS 
339 %This works with one frequency & one omega 
340 if choice2== 
341 if choice3== 
342 if choice4==1 %Fixed Free 
343 if modenumber== 
344 mode=input (‘input mode number you are checking\n'); 
345 end 
346 DD(:,mode)=theta(:,1,1); 
347 GG(:,mode)=T(:,1,1); 
348 sign=(-1)%* (mode+l1); 
g 





[rowgage, gages]=size(span) ; 

Eliminate a Bad Strain Gage??? 

351 if eliminate== 

352 eliminate=input ('enter 1 if you want to zero out a strain gage or 
zero not to eliminate any\n'); 

353 if eliminate== 

354 maxmodes=gages-2; 





Ww 
ol 
jo) 
ole 














355 disp(['If the number of modes calculated exceed ', 
num2str(maxmodes) ,' your result will be incorrect']); 
356 killit=input('enter the strain gage location number\n (closest to 
the root = 1) you want to eliminate\n'); 

357 for g=l:gages 

358 if g<killit 

359 span (rowgage,g)=span (rowgage,gQ) ; 

360 end 

361 aif g>killit 

362 span(rowgage,g-1)=span (rowgage,g) ; 

363 end 

364 end 

365 gages=gages-1; 

366 span=span (rowgage, 1l:gages) 

367 end 

368 end 

369 %%SS5SSSEND OF ELIMINATION ROUTINESS%3%%%S%S%% 











370 %SInterpolate calculated harmonic torques to span locations for 
direct comparison 
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371 





for kkk=nodes:-1:2 






























































(closest to 


372 for ccc=l:gages 

373 if x(kkk)>=span(ccc) 

374 twotorque=[GG(kkk,mode) GG(kkk-1,mode) ]; 

375 twox=[x(kkk) x(kkk-1)]; 

376 coeff=polyfit (twox,twotorque, 1); 

377 G(ccc,mode) =coeff(1).*span(ccc) +tcoeff (2); 

378 end 

379 ccc=ccct1; 

380 end 

381 end 

382 xl=0:1; 

383 yl=zeros(size(xl)); 

384 figure 

385 if choice== 

386 plot(x,theta(:,1,1),'ks',x, (Sign*sin((mode-0.5) *pi*x)),'k-') 
387 end 

388 if choice~=2 

389 if choice~=5 

390 plot(x,theta(:,1,1),'k-s"') 

391 end 

392 end 

393 if choice== 

394 plot(x(3:nodes),theta(3:nodes,1,1),'k-s"') 

395 end 

396 line(xl,yl) 

397 xlabel('r/R') 

398 ylabel('Torsional Deformation (rad)') 

399 title(['Mode Shape # ',num2str(mode),'-- ',num2str (nodes), ' 
Segments -- ','Wn = ',num2str(w),' rad/sec']) 

400 legend('calculated mode shape', 0) 

401 end %Sfor choice4==1 

402 if choice4==0 %FREE FREE 

403 if modenumber== 

404 mode=input (‘input mode number you are checking\n'); 
405 end 

406 DD(:,mode)=theta(:,1,1); 

407 GG(:,mode)=T(:,1,1); 

408 sign=(-1)%* (mode+l1); 

409 [rowgage, gages]=size(span) ; 

410 %Eliminate a Bad Strain Gage??? 

411 if eliminate== 

412 eliminate=input ('enter 1 if you want to zero out a strain gage or 
zero not to eliminate any\n'); 

413 if eliminate== 

414 maxmodes=gages-2; 

415 disp(['If the number of modes calculated exceed 
',num2str(maxmodes),' your result will be incorrect']); 
416 killit=input('enter the strain gage location number\n 
the root = 1) you want to eliminate\n'); 

417 for g=l:gages 

418 if g<killit 

419 span (rowgage,g)=span (rowgage,gQ) ; 

420 end 

ADL af g>killit 

422 span(rowgage, g-1)=span(rowgage,g) ; 

423 end 

424 end 

425 gages=gages-1; 

426 span=span (rowgage, l:gages) 

427 end 

428 end 

429 %%%%S%SSSEND OF ELIMINATION ROUTINES%S%%%%%%3% 
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430 
direc 
431 
432 
433 





SInterpolate calculated harmonic torques to span locations for 
t comparison 

if choice~=3 

for kkk=nodes:-1:1 

for ccc=l:gages 

if x(kkk) >=span (ccc) 

twotorque=[GG(kkk,mode) GG(kkk-1,mode) ]; 

twox=[x(kkk) x(kkk-1)]; 

coeff=polyfit (twox,twotorque,1); 

G(ccc,mode) =coeff(1).*span(ccc)+coeff (2); 











end 

ecc=ccctl; 

end 

if kkk>1 

kkk=kkk-1; 

end 

end 

end 

if choice~=2 

plot (x,theta(:,1,1),'k-s') 

end 

if choice==2 %plots DNW from offset to tip 
plot (x(3:nodes),theta(3:nodes,1,1),'k-s"') 
end 

line (xl,yl) 

xlabel('r/R"') 

title(['Mode Shape #',num2str (mode), '-—',num2str(nodes),' Segments 
"Wn = ',num2str(w),' rad/sec']) 

if choice~=3 

ylabel('Torsional Deformation (rad)') 
legend('calculated mode shape',0) 

end 

end Sfor choice4==0 loop 

end %for choice3==1 line 

end %for choice2==1 line 





if modenumber== 

ccontinue=input ('enter one to try an exact frequency to check or 
to quit\n'); 

end 

if ccontinue== 

if modenumber== 

w=input ('enter frequency to check in rad/sec\n'); 

omega=input ('enter rotational speed to check in rad/sec\n'); 

end 

choice3=1; 

choice2=1; 

upper=0; 

lower=0; 

upperomega=0; 

loweromega=0; 

end 

if modenumber== 

if ccontinue~=1 

ccontinue=input ('press zero if you are sure you want to quit, else 
1?') 





modenumber~=0 
mode==modenumber 
continue=0; 


Fh Fh 





f mode<modenumber 
mode=modet+1; 
w=rootfreq (mode) ; 
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490 
491 
492 
493 
494 
495 
496 
497 
498 
499 


choice3=1; 
choice2=1; 
upper=0; 
lower=0; 
upperomega=0; 





loweromega=0; 
ccontinue=1; 
end 

end 

end 


500 modenumber=mode; 
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APPENDIX B. STRAIN PATTERN ANALYSIS PROGRAM 


A. INSTRUCTIONS FOR RUNNING THE SPA PROGRAM 


The option to run the SPA program is given after the Holzer program generates 
the number of mode shapes desired. Your initial choice of how many mode shapes to 
generate cannot be changed unless you start over and run the Holzer program again. If 
you choose to run the SPA program line 4 gives the option to see plots of SPA estimated 
torsion moments versus measured torsion moments for each harmonic (lines 68 to 78). 
Similarly, line 5 of the program gives the option to plot the SPA vector by mode number 
(lines 81 to 154). These options are not recommended if you are only concerned with the 
final estimated deflection result of the SPA program. You will also have the option to 
choose to estimate the first 10 harmonics (21 terms through the 10th sine) or 20 
harmonics (41 terms through the 20th sine). The matrix of measured torsion moments 
you enter must be equal to or greater than your choice of harmonics to estimate. Enter 
the measured torsion moments in the same units that correspond to those you chose in 
Table 15. The columns of your measured torsion moment matrix are the harmonics and 
the rows are the measurement stations. After entering the measured moments you can 
opt to randomly vary the measured data to replicate measurement error. Line 20 asks for 
a vector of span locations of comparison deflection data to compare against the SPA 
estimated deflection. Lines 24 to 37 eliminate an entire set of strain gage data from the 
measured moment matrix if you selected this option during the Holzer program. Lines 39 
to 48 interpolate the calculated mode shapes to the same span locations of the comparison 
deflection values entered in line 20. Line 62 is Step One of the SPA method using 
Equation 1.7 to calculate the SPA vector. Line 63 is Step Two of the SPA method using 
Equation 1.8 to calculate the estimated torsion moment matrix to compare against the 
measured torsion moment matrix. Line 64 is Step Three of the SPA method using 
Equation 1.9 to calculate the estimated deflection. All three Equations are inside a loop 
that runs once for each harmonic term. Each SPA vector is stored as a column in a 
matrix of {q}. Similarly, each harmonic of the estimated deflection solution {deflsol} 
and estimated moments {Checkq} are stored as columns of their respective matrix. The 


remainder of the program produces plots based on the earlier choices and converts the 
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harmonics of measured torsion moments and estimated deflections (radians) into the time 
domain. 


B. SPA PROGRAM VARIABLES 


steadytendefl The (nx1024) matrix of estimated SPA deflection (radians) where each 
row corresponds to the span locations entered for the comparison data and the 1,024 


columns range from zero to 27 radians (line 193). 


deflsol The (nxh) matrix of harmonics of SPA estimated deflection where the 
rows correspond to the interpolation stations chosen for comparison and the columns 


correspond to harmonic terms (line 64). 


newxaxis The final SPA estimated deflection in deflsol or steadytendefl may not be 
available at all comparison stations if the last calculated mode shape location did not 
exceed the last comparison deflection location. This variable contains all of the span 
locations that the SPA deflection values can be compared to the chosen set of deflections 


(line 67). 


steadyten The (gx1024) matrix of measured torsion moment values (ft-lb) where 
each row corresponds to the measurement span locations and the 1,024 columns range 


from zero to 2 7 radians (line 173). 


M The (gxh) matrix of harmonic terms for the measured torsion moment 


entered by the user (ft-lb or N-m) (line 16). 


xaxis The user input span locations to interpolate values of the torsional 


deflection mode shape [DD] from the Holzer program (non-dimensional, r/R) (line 20). 


D The (xaxis,m) matrix of Hozler calculated deflection mode shapes 
interpolated to the desired span locations by row and each column corresponds to a 


separate mode (line 47). 


q The (mxh) matrix of SPA vectors where each column is the SPA vector 


for each harmonic, and the rows correspond to mode numbers (line 62). 


checkq The (gxh) matrix of SPA estimated torsion moments where rows 


correspond to measurement span locations and columns correspond to harmonic terms 
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GC; SPA PROGRAM LISTING 


1 SSTART OF SPA PROGRAM, PART 2% 

2 GOSPA=input (‘enter 1 if you want to conduct Strain Pattern 
Analysis\n'); 

3 if GOSPA== 

4 figyes=input('press 1 if you want to see a plot of Test Data vs SPA 
Torsion Moment vs. r/R by Harmonic\n'); 

5 figyes2=input ('press 1 if you want to see the proportions (q) SPA 
vector plotted vs. mode number by harmonic\n'); 

6 harmonic=input (‘enter the 10 (through 10th Sine) or 20 (through 
20th Sine) for the number of harmonics you want to solve\n'); 

7 harmonicterm=2*harmonictl; 

8 if harmonic==10 

9 fortyone=0; 

10 end 

11 if harmonic==20 

12 fortyone=1; 

13 end 

14 if choice~=0 

15 if choice~=2 

16 M=input ('enter the harmonic torsion moment (FT-LBS)\n'); 

18 STHE NEXT inputs ARE REQUIRED if you want to plot you're results 
































against other results?? 








19 SLocations entered below are used to interpolate deflection mode 
shape data for final SPA result to be easily compared to other methods 
20 xaxis=input ('enter the vector of r/R locations for the comparison 
set of deflection information\n'); 

Zi end 

22 end 

23 rowM colM]=size(M); 

24 SIf you chose to eliminate a Bad Strain Gage??? 

25 if eliminate== 

26 for g=1:rowM 








27 if g<killit 
28 M(g,:)=M(g,:); 





























29 end 

30 aif g>killit 

31 M(g-1, )=M(g, )e 

32 end 

33 end 

34 rowM=rowM-1; 

35 M=M(1:rowM,:); 

36 end 

37 S$SSSSSEND OF ELIMINATION ROUTINES%%3%%%%%% 
38 [rowxaxis colxaxis]=size(xaxis); 

39 SInterpolate mode shape to xaxis locations to be able to compare 
test to SPA 

40 for mmode=1:mode 

4l for kkk=Nodes:-1:1 

42 for ccc=1:colxaxis 

43 if x(kkk)>=xaxis (ccc) 


44 twodefl=[DD(kkk,mmode) DD(kkk-1,mmode) ]; 
45 twox=[x (kkk) x(kkk-1)]; 
46 coeff=polyfit (twox,twodefl,1); 


47 D(ccc,mmode) =coeff(1).*xaxis(ccc)+tcoeff (2); 
48 end 

49 ecc=ccctl1; 

50 end 


51 if kkk>1 
52 kkk=kkk-1 
53 end 
54 end 
55 mmode=mmode+1; 





~“e 
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56 end 


Deh SSize of G (Referred to as [T] in thesis) is (#r/R stations of 
Strain Gages x # modes in solution (#)) 

58 Size of Mis (#r/R Stations in Test Data x # Harmonic Terms (21 
or 41)) 

59 Size of q is #modes x # Harmonic terms 

60 for amplitudes=1:harmonicterm; 


61 STHIS IS THE SPA Vector {q} with number of modes used for number of 








62 q(amplitudes) =inv(G.'*G)*G.'*M(:,amplitudes) ;%Equation 1.7 Step One 
63 checkq(:,amplitudes) =G*q(:,amplitudes) ;%Equation 1.8 Step Two 








64 deflsol(:,amplitudes) =D*proportions(:,amplitudes);%Equation 1.9 
Step Three 

65 [rowD colD]=size(D); 

66 nodes=rowD; 

67 newxaxis=xaxis (1l:nodes) ; 


68 if figyes== 

69 x1=0:1; 

70 yl=zeros(size(xl)); 

71 figure 

72 plot (span, checkq(:,amplitudes),'k-d',span,M(:,amplitudes),'b:+') 
73 line (xl,yl) 

74 Stitle([num2str(modenumber) ,' mode solution: Test Data vs SPA 
Torsion Moment for Harmonic Term ',num2str (amplitudes) ]) 

75 xlabel('r/R') 

76 ylabel('Torsion Moment (ft-lbs) ') 












































77 Slegend('SPA Torsion Moment', 'Measured Torsion Moment', 0) 
78 end 

719 amplitudes=amplitudest1; 

80 end 

81 if figyes2== 

82 xl=1:mode; 

83 yl=zeros(size(xl)); 

84 figure 

85 plot (1:1:mode,q(:,1),'k-d',1:l1:mode,q(:,2),'b- 
o',1:1l:mode,q(:,3),'g-*') 

86 line (xl,yl) 

87 xlabel('mode number') 

88 ylabel('SPA Vector") 

89 title([num2str(modenumber) ,' mode solution: Steady, lst Cos, and 
lst Sin Least Squares Proportions']) 

90 legend('Steady','lst Cos','1lst Sin',0O) 

91 figure 

92 plot (1:1:mode,q(:,4),'k-d',1:l:mode,q(:,5),'b-o') 

93 line (xl,yl); 

94 xlabel('mode number') 

95 ylabel('SPA Vector") 

96 title ([num2str(modenumber) ,' mode solution: 2nd Cos and 2nd Sin 
Least Squares Proportions']) 

97 legend('2nd Cos','2nd Sin',0) 

98 figure 


99 plot (1:1:mode,q(:,6),'k-d',1:l:mode,q(:,7),'b-o') 

100 line(xl,yl) 

101 xlabel('mode number') 

102 ylabel('SPA Vector') 

103 title([num2str(modenumber) ,' mode solution: 3rd Cos and 3rd Sin 
Least Squares Proportions']) 

104 legend('3rd Cos','3rd Sin',0O) 

105 figure 

106 plot(1:l:mode,q(:,8),'k-d',1:1:mode,q(:,9),'b-o') 

107 line(xl,yl) 

108 xlabel('mode number') 

109 ylabel('SPA Vector") 
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10 title([num2str(modenumber) ,' mode solution: 4th Cos and 4th Sin 
east Squares Proportions']) 
legend('4th Cos','4th Sin',0) 
figure 
plot (1:1:mode,q(:,10),'k-d',1:l:mode,q(:,11),'b-o') 
line (xl,yl) 
xlabel('mode number') 
ylabel('SPA Vector') 
if title ([num2str(modenumber) ,' mode solution: 5th Cos and 5th Sin 
Least Squares Proportions']) 
118 legend('5th Cos','5th Sin',0O) 
119 figure 
120 plot (1l:l:mode,q(:,12),'k-d',1:1:mode,q(:,13),'b-o') 
121 line(xl,yl) 
122 xlabel('mode number') 
123 ylabel('SPA Vector") 
124 title([num2str(modenumber) ,' mode solution: 6th Cos and 6th Sin 
Least Squares Proportions']) 
125 legend('6th Cos','6th Sin',0O) 
126 figure 
127 plot (1l:1l:mode,q(:,14),'k-d',1:1:mode,q(:,15),'b-o') 
128 line(xl,yl) 
129 xlabel('mode number') 
130 ylabel('SPA Vector') 
131 title([num2str(modenumber) ,' mode solution: 7th Cos and 7th Sin 
Least Squares Proportions']) 
132 legend('7th Cos','7th Sin',0) 
133 figure 
134 plot(1l:l:mode,q(:,16),'k-d',1:1:mode,q(:,17),'b-o') 
135 line(xl,yl) 
136 xlabel('mode number') 
137 ylabel('SPA Vector") 
138 title([num2str(modenumber) ,' mode solution: 8th Cos and 8th Sin 
Least Squares Proportions']) 
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139 legend('8th Cos','8th Sin',0O) 

140 figure 

141 plot(1l:l:mode,q(:,18),'k-d',1:1:mode,q(:,19),'b-o') 
142 line(xl,yl) 

143 xlabel('mode number') 

144 ylabel('SPA Vector") 

145 title([num2str(modenumber) ,' mode solution: 9th Cos and 9th Sin 
Least Squares Proportions']) 

146 legend('9th Cos','9th Sin',0O) 

147 figure 

148 plot(1l:l:mode,q(:,20),'k-d',1:1:mode,q(:,21),'b-o') 
149 line(xl,yl) 


150 xlabel('mode number') 

151 ylabel('SPA Vector") 

152 title([num2str(modenumber) ,' mode solution: 10th Cos and 10th Sin 
Least Squares Proportions']) 

153 legend('10th Cos','10th Sin',0O) 

154 end 

155 %$%%SConvert from Torsional Moment Harmonics to Time Domain%%% 
156 sta=1; 

157 for sta=1:rowM 

158 MM=M(sta,:); 

159 Harmtors=MM(:); 

160 harm=1; 

161 for harm=1:harmonic; 














162 A(harm) =—-1*Harmtors (2*harm) ; 
163 B(harm) =-1*Harmtors (2*harmt1) ; 
164 i=1; 


165 for i=1:1024 
166 Lattors (harm, i)=A(harm) *cos (harm*CNEWAZIMUTH (i) ); 
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167 Longtors (harm, i)=B (harm) *sin (harm*CNEWAZIMUTH (1) ); 
168 end 

169 SNow put collective and cyclic together Harmonic 
170 Tors (harm, :)=-Lattors (harm, :)-—-Longtors (harm,:); 
171 harm=harm+1; 

172 end 

173 steadyten(sta,:)=M(sta,1)+sum(Tors,1); 

174 sta=statl; 

175 end 








176 S%SNOW Convert DEFL Harmonics to TIME HIST%%% 
177 for sta=1:rowD 

178 ff=deflsol(sta,:); 

179 Harmdefl=ff(:); 

180 harm=1; 

181 for harm=1:harmonic; 

182 A(harm) =-1*Harmdefl (2*harm) ; 

183 B(harm) =—-1*Harmdefl (2*harmt+1) ; 

184 i=1; 




















185 for 1i=1:1024 

186 Latdefl (harm, i)=A(harm) *cos (harm*CNEWAZIMUTH (i) ); 
187 Longdefl (harm, i)=B (harm) *sin (harm*CNEWAZIMUTH (1) ); 
188 end 

189 SNow put collective and cyclic together Harmonic 
190 DEFL(harm, :)=-Latdefl (harm, :)-Longdefl(harm,:); 
191 harm=harmt1; 

192 end 

193 steadytendefl(sta,:)=deflsol(sta,1)+sum(DEFL,1);%Radians 
194 sta=statl; 

195 end 

196 end %end for GOSPA== 
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APPENDIX C. UTRC MODEL-SCALE BLADE PROPERTIES 


0.18600 
0.21601 
0.24000 
0.39100 
0.42600 
0.46600 
0.49801 
0.55901 
0.59201 
0.66601 
0.68301 
0.73601 
0.78200 
0.78601 
0.81301 
0.85800 
0.86200 
0.86500 
0.88260 
0.90260 
0.91300 
0.92891 
0.93250 
0.93610 
0.95391 
0.95611 
0.97390 
0.98110 


0.23401 
0.27801 
0.40901 
0.46401 
0.49800 
0.55900 
0.59200 
0.66600 
0.68300 
0.73600 
0.76601 
0.78600 
0.81300 
0.85400 
0.85931 
0.86221 
0.88001 
0.90100 
0.91111 
0.92890 
0.93201 
0.93601 
0.95390 
0.95610 
0.97200 
0.97401 





rossi | 1.00000, | 
Table 16. UTRC Blade 126-Segment Span Locations. [After: Ref. 17] 
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152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 152.78 
152.78 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 72.22 
72.22 

ea ee a 

Table 17. .UTRC Blade Torsional Stiffness. [After: Ref. 17] 
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1.79E-05 
8.29E-05 1.79E-05 
1.79E-05 1.86E-05 
1.86E-05 6.22E-05 
6.22E-05 6.22E-05 
3.63E-05 6.48E-05 
6.48E-05 6.48E-05 
6.48E-05 4.14E-05 
7.25E-05 7.25E-05 
7.25E-05 7.25E-05 
5.61E-05 5.96E-05 
4.14E-05 4.14E-05 
4.14E-05 2.11E-05 
2.11E-05 7.10E-05 
7.10E-05 3.11E-05 
3.11E-05 5.18E-05 
5.18E-05 1.04E-04 
1.04E-04 1.04E-04 
3.11E-05 3.09E-05 
3.08E-05 1.76E-05 
1.76E-05 9.98E-05 
0.00E+00 0.00E+00 
1.44E-04 1.08E-04 
1.02E-04 
fO.00E+oo]0.o0E+ooT—CdTCCd 
Table 18. .UTRC Blade Torsional Mass Moment of Inertia. [After: Ref. 17] 
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5.93E-11 5.86E-11 
1.14E-10 6.68E-11 
6.66E-11 6.52E-11 
5.99E-11 6.65E-11 
6.54E-11 5.25E-11 
5.24E-11 5.03E-11 
5.03E-11 8.35E-11 
5.15E-11 4.87E-11 
4.87E-11 4.17E-11 
6.53E-11 6.64E-11 
6.64E-11 2.98E-11 
2.09E-11 2.80E-11 
2.63E-11 2.14E-11 
2.14E-11 1.90E-11 
1.90E-11 8.77E-12 
8.55E-12 6.19E-12 
6.08E-12 6.07E-12 
5.96E-12 8.97E-12 
8.78E-12 7 A8E-12 
7.36E-12 6.75E-12 
6.61E-12 5 45E-12 
5.45E-12 2.06E-12 
1.67E-12 2.30E-12 
2.30E-12 1.68E-12 
1.68E-12 1.64E-12 
1.64E-12 9.73E-13 


9.08E-13 | 9.08E-13 | 9.05E-13 | 4.45E-13 
3.25E-13 | 3.25E-13 | 1.55E-13 | 1.91E-14 
1.89E-14 | 0.00E+00 


Table 19. _ UTRC Blade Tension-Torsion Term without Square of Rotational Speed. [After: 
Ref. 17] 
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| 
P9325 9.00] 8.7] 8.87 
[6.83] 6.96] 6.00] 7.37 
| 8.82] BOP 

Table 20. UTRC Blade Feathering Angle Before Deformation. [After: Ref. 14] 
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APPENDIX D. STRAIN PATTERN ANALYSIS PLOTS 


UH-60A Model Rotor, DNW Run 13.20, 4“ =0.301 


— Steady 


-&- 1st Cos 
—#— 1st Sin 


SPA Yector 





4 2 3 4 5 
mode number 


Five-Mode SPA Vector for the Steady, 1st Cosine and Ist Sine Terms. 





Figure 44. 
é 
> 
at 
a) 
= 2 3 4 5 
mode number 
Figure 45. Five-Mode SPA Vector for the 2nd Cosine and 2nd Sine Terms. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 


SPA Yector 


—@ 3rd Cos 
-@- 3rd Sin 





“4 2 3 4 5 
mode number 


Figure 46. Five-Mode SPA Vector for the 3rd Cosine and 3rd Sine Terms. 


—6 4th Cos 
—& 4thSin 


0.015 





SPA, Vector 


2 3 4 5 
mode number 


Figure 47. Five-Mode SPA Vector for the 4th Cosine and 4th Sine Terms. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 


— 5th Cos 
-@- Sth Sin 


SPA, Yector 





“4 2 3 4 5 
mode number 


Figure 48. Five-Mode SPA Vector for the 5th Cosine and 5th Sine Terms. 


SPA, Yector 





“4 2 3 4 5 
mode number 


Figure 49. Five-Mode SPA Vector for the 6th Cosine and 6th Sine Terms. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 


SPA Yector 


—6- 7th Cos 
= 7thsSin 





“4 2 3 4 5 
mode number 


Figure 50. Five-Mode SPA Vector for the 7th Cosine and 7th Sine Terms. 
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Figure 51. Five-Mode SPA Vector for the 8th Cosine and 8th Sine Terms. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 52. Five-Mode SPA Vector for the 9th Cosine and 9th Sine Terms. 


SPA Yector 


—6- 10th Cos 
-= 10th Sin 


“4 2 3 4 5 
mode number 





Figure 53. Five-Mode SPA Vector for the 10th Cosine and 10th Sine Terms. 
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Figure 54. 
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UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 
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UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 
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Figure 56. One-Mode SPA versus Measured Torsion Moments for 8th Sine to 10th Sine. 


102 


Figure 57. 
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Figure 58. | Two-Mode SPA versus Measured Torsion Moments for 4th Sine to 8th Cosine. 
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Figure 59. 
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Figure 60. 


UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 
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Three-Mode SPA versus Measured Torsion Moments for Steady to 4th Cosine. 


106 


14 





= 
& 


Torsion Moment (ft-lbs) 
° 
& 


° 
= 


O2F 





1st COS 








O41 


O41 


08 og 





2nd COS 
+ 





-041F 


Torsion Moment (ft-lbs) 





0.25 
0 





O41 


08 og 





=] ° 
a is 


Torsion Moment (ft-lbs) 
° 
ES 


-06F 





3rd COS 








O41 


02 


03 


o4 


Os 
oR 


O68 


O7 


os 


09 





O04 


02 


0.3, 


-04 


Torsion Moment (ft-lbs) 


05 


06 








4th COS 


+ 4 








07 
0 


o4 


04 


0s 
oR 











UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 









































































































































as a4 
+ 
4th SIN . Sth COS 
ost 035 : 4 
ost 4 
oat 
z Bo2s 4 
= 03 & 
5 5 
e & 02- 4 
5 5 
= = 
§ 02r § 
g So15- 4 
& & 
ott 
aap 4 
9 * ost 4 
04 1 , 1 . . 1 : A . 0 . 1 1 1 1 1 fl 1 1 
of o2 03 o4 o5 O8 o7 o8 o9 1 o of 02 O38 O4 OS O08 OF oO8 os 1 
MR oR 
02 o04 
Sth SIN 
oast 
oo2b 4 
otf if 
= a 
3B a 
8 8 
= 00st = 
3 a. 3 
2 2 -0.02+ 4 
5 5 
= = 
5 9 5 
3 3 
5 5 
4 © 004+ 4 
-0.05+ 
oak -0.08+ 4 
“048 A 1 1 1 1 -0.08 1 1. 1 A A 1. A +) 
ood os o7 o8 a9 1 o of 02 O03 O4 OS O08 OF oO8 og 1 
oR 
002 1 r : 0.02 r r : r : 
+ 
oof th SIN : 4 7th COS 
+ 
a 
a 
-ootf 4 
a 0.02 4 
a a 
= 002} 4 = 
5 5 
& -0.03+ 4 & -0.04¢ 4 
5 5 
= = 
F F 
$ -0.04- 4 5 
e ® oost | 
-0.05 4 
-0.08 4 
-0.08+ 4 
-o07b + 4 
-0.08 04 * 
o of 02 O03 O4 OS O08 o7 o8 og 1 o of o2 O03 Of OS 08 o7 o8 o9 1 
oR oR 
0.05 0.08 
7th SIN Gace 
+ 
7 0.06 + 4 
-0.05 4 
0.04 4 
B o4 ] 3 
= = 002 4 
a a 
€ -0.15 4 13 
5 5 
= = 
5 5 ° 
3 3 
eee 1 & 5 
-0.02 4 
0.25 4 
oS 
“n3 a 0.04 4 
035 1 +: . . 1 . . 1 0.06 1 . A 1 A 1 1 . . 
oof 02 03 of 05 O08 OF O08 os 1 oof o2 O03 Of OS O08 OF o8 o9 1 
oR oR 


Figure 61. Three-Mode SPA versus Measured Torsion Moments for 4th Sine to 8th Cosine. 
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UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 
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Figure 62. Three-Mode SPA versus Measured Torsion Moments for 8th Sine to 10th Sine. 


108 


Figure 63. 
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UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 
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Figure 64. Four-Mode SPA versus Measured Torsion Moment for 4th Sine to 8th Cosine. 
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Figure 65. 
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UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 
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Figure 66. Five-Mode SPA versus Measured Torsion Moments for Steady to 4th Cosine. 
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UH-60A Model Rotor, Measured = Dotted Line, SPA = Solid Line 
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Figure 67. Five-Mode SPA versus Measured Torsion Moments for 4th Sine to 8th Cosine. 
113 


Figure 68. 
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Five-Mode SPA versus Measured Torsion Moments for 8th Sine to 10th Sine. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 69. Five-Mode SPA versus Measured Torsion Moment Time History at x = 0.2. 


30 





Torsional Moment (in-Ib) 


o 50 100 150 200 250 300 350 
Azimuth (deg) 


Figure 70. Five-Mode SPA versus Measured Torsion Moment Time History at x = 0.35. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 71. Five-Mode SPA versus Measured Torsion Moment Time History at x = 0.5. 
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Figure 72. Five-Mode SPA versus Measured Torsion Moment Time History at x = 0.642. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 73. Five-Mode SPA versus Measured Torsion Moment Time History at x = 0.745. 
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Figure 74. Five-Mode SPA versus Measured Torsion Moment Time History at x = 0.875. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 75. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.225. 
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Figure 76. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.4. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 77. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.55. 
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Figure 78. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.675. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 79. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.775. 
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Figure 80. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.865. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 81. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.92. 
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Figure 82. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.945. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 83. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.965. 
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Figure 84. Five-Mode SPA versus DNW Calculated Deformation Time History at x = 0.99. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 85. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.225. 


15 


So 
in 


SPA Deflection (deg) 
o 
wn o 


-1 * steady 





-2 15 -1 -05 0 os 1 15 2 
DNVY Deflection (deg), b= -0.021183, corr. coef. = 0.95916, slope= 1.2937 


Figure 86. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.4. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 87. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.55. 
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Figure 88. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.675. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 89. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.775. 
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Figure 90. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.865. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 91. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.92. 
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Figure 92. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.945. 
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UH-60A Model Rotor, DNW Run 13.20, “=0.301 
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Figure 93. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.965. 
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Figure 94. Linear Regression Plot of Five-Mode SPA versus Measured Torsional 
Deformation by Harmonic at x = 0.99. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 95. Linear Regression Y-intercept versus Span. 
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Figure 96. Linear Regression Correlation Coefficient versus Span. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 97. Linear Regression Slope versus Span. 
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Figure 98. Varied Structural Property versus Original Five-Mode SPA at x=0.400. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 








Torsional Deflection (deg) 














Original SPA 
—— Varied SPA 
i 1 i i 1 1 
0 50 100 150 200 250 300 350 
Azimuth (deg) 





Figure 99. Varied Structural Property versus Original Five-Mode SPA at x=0.775. 
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Figure 100. Eliminated Strain Gage Four-Mode versus Original Five-Mode SPA at x=0.400. 
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UH-60A Model Rotor, DNW Run 13.20, “~=0.301 
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Figure 101. Eliminated Strain Gage Four-Mode versus Original Five-Mode SPA at x=0.775. 
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Figure 102. Varied Measurement versus Original Five-Mode SPA at x=0.400. 
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UH-60A Model Rotor, DNW Run 13.20, “4 =0.301 
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Figure 103. Varied Measurement versus Original Five-Mode SPA at x=0.775. 
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